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ON DIOPHANTINE APPROXIMATION AND THE 
GENERALIZED ELLIPTIC THETA FUNCTION 


H. 8. A. POTTER (Aberdeen) 
[Received 1 November 1937] 


1. Introduction 
m m 
Let Q(x) => }> 8,2;2x, be a positive definite quadratic form with 


j=1 k=1 
integer coefficients and 


f(r) = Der" (r = ax+Iy; y > 0), (1) 


where the summation is over all integer values of x,,..., z,,. It is well 
known that the elliptic theta function f(r) defined by this series is 
analytic throughout the half-plane y > 0 and that the line y = 0 is 
a natural boundary of the function. We shall investigate the 
asymptotic behaviour of f(r) as 7 moves up to various types of points 
on the natural boundary. 

Let A(n) denote the number of integer solutions of 


Q(x) =n. 
Then, as the series (1) is absolutely convergent for y > 0, we may 
write it in the form 


f(r) = 1+ > A(nje7i*™, (2) 
n=1 
This is the Abel sum of the trigonometric sum 


N 

R(N,x) = > A(nje7™ 
n=0 

—_ > green, (3) 

Q(E)SN 
whose asymptotic behaviour we shall also consider. 

The solution by Hardy and Littlewood* of the corresponding 
problem for © 
3,(0|r) = adel 

n=-—® 
depends on the transformation formulat 
ar+b 
cr+d 


3,(0|7) = «(er+d)- 'a,(0 


(«x a constant), 


* Hardy and Littlewood, 1. + See Hardy and Littlewood, 1, p. 227. 
3695.9 M 

















162 H. 8. A. POTTER 
which holds for transformations of the modular group. Siegel* has 


proved a similar formula for f(r) but only for the sub-group of the 
modular group satisfying the congruence 


a b 1 0 _ 
= 28), 
(’ i} (; ) (mod 28) 


~ 


where S denotes the determinant of the matrix S = (s,,). This con- 
dition makes it useless for our purpose. We shall prove a formula 
(see Lemma 4) less precise than Siegel’s, but sufficient for our purpose 
and free from such conditions. Otherwise, the proof follows the 
method of Hardy and Littlewood. 


2. Notation 

We shall adopt the notation of Siegel. If & is a matrix, let Y’ 
denote the transpose of Y, ie. the matrix obtained from Y by inter- 
changing rows and columns. Thus G’ = GS for the symmetric matrix 
S. In the case when % is a non-singular square matrix we denote its 
inverse by U-!. The column of elements 2, 2g,..., 2, will be denoted 
by the corresponding small German letter x, so that x’ will be the 
row of elements (%,,...,2,,). We shall call a matrix an integer when all 
its elements are integers. If U is an integer matrix whose determinant 
|U| = 1, we shall call it a unit matrix. We shall write 

2 = B (mod q) 

when all the elements of U—%B are divisible by q. 

To indicate a summation over all integer columns x we shall write 


>, while >’ will mean such a summation with the null column n 
¥ 


¥ 
omitted. For a summation over a complete remainder system to 


modulus 6 we shall write > (i.e. each element of x runs over a com- 
x(b) 


plete remainder system to modulus b). If x runs over all those integer 
columns of a complete remainder system to modulus b which satisfy 
the congruence %z =n (mod), where 8|b, we shall write > 


x(b) 
Ax=n(P) 
We reserve S to denote the positive definite symmetric integer 
matrix (s,;,) of m rows and columns and determinant S = |S|. Then 


the quadratic form under discussion may be written 
Q(x) = x’ Gx. 
If all the elements of the principal diagonal of S are even, then Q(x) 


* Siegel, 1, Lemma 30. 
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can only represent even integers. We therefore describe S as even in 
this case; otherwise we say that © is odd. 

Let d,, be the greatest common divisor of all the sub-determinants 


of S of the kth order, so that d,, = |S|. Write 
é, = d,, e, =d,/d,_, (k = 2,3,...,m) 
Then the numbers ¢é,, ég,..., €@,, Which may easily be shown to be 


integers, are called the elementary divisors of G. 


3. Congruences 
The central theorem in the theory* of elementary divisors is 
LeMMA 1. There exist two unit matrices U and B such that 
yu-'SB = D, 
where D is the diagonal matrix (e;8;;), 5;; being Kronecker’s symbol. 
An easy deduction from Lemma 1 is 


Lemma 2.} Let N(S,b)= D> 1, 
x(b) 
Sr=n(b) 
i.e. N(S,b) is the number of integer solutions of 
Sx = n (modb) 
which are incongruent to modulus b. Then 


N(G,b) = 8, 8...8, 


where 8), = (€;, 6). 
We note that, as s, < é,, 
N(G, 6) < €,¢...4, = d,, = B. (4) 
Also, if 6 is odd, the congruence 2Sx = n (modb) implies Sz = n 
(modb)sothat wos 5) = N(G,b) (b odd). (5) 
If b is even and b = 28, the congruence 2Szx = n (modb) implies 
Sx = n(mod 8) so that, using the substitution x = By+3, 


N26,b)= FY l= SY Y1=2"N(G,A). (6) 
_ x(b) _ 3(B) . 92) 
Sr=n(f) S3=n(f) 


In either case we have 
N(26,6) < 2"8. (7) 


* See H. J. S. Smith 1, 391, or Bachmann, 1. 
+ See e.g. H. J. 8. Smith, 1, 399. 
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4. The generalized Gaussian sum 


Let us now consider the generalized Gaussian sum 


G(S, a/b, bh) = > e27ialby’Gr+4'}, 


x(b) 
where (a,b) = 1, b > 0, and § is a column of m rational elements such 
that bh is an integer column. When § = n, the null column, we write 

G(S, a/b, n) = G(S, a/b). 
The following lemma is merely a refinement on that of Siegel.* 


Lemma 3. We have 
G(S, a/b, h)| < b™#{N(26, b)}* < 2mi2Sipm?, 


and 
[ bm2{N(S, b)}! b odd), 


( 
0 (b = 28; B odd and & odd), 
|G(S, a/b) | (2by"2{N(S,B)}! (b = 2B; B odd and S even), 
0 ( 
(2b)m2EN(S, B)} 


5 
b = 28; B even and € odd), 
b = 2B; B even and € even), 


| 
| 
| 
l 


where 
€ = (cy) = (Es vn); W = (wy) = VU, o;, = (e;; 8). 
o; O;; 
The matrix € is clearly an integer. By Lemma | 
SB =uUD 

so that WD = PUD = VSB = (VB’SB)’ = (WD), 

e. WD (= (w;,,e;,)) is symmetric. Hence € is symmetric. Let 2° 
2px, 2k, 2% be the highest powers of 2 dividing f, e;, o;,, W,);, respec- 
tively. Then the highest power of 2 which divides c,, is 27*, where 


Ve = |P—PR| TOK (8) 
Thus € is even if |p—p,|+o; > 1 for k = 1, 2,..., m, and € is odd if 
Py = p and a, = 0 for some value of k. 


Proof of Lemma 3. Writing r = a/b, we have 


G(S, 7; b)|? = G(S, 7, 5)G(S, —f, —h) 
= b >on Ti{r(3’S3—y’Sy)-+(3'—y’)b} 
a(b) v(b) 


on Oe e27i{re’Sx+2ry’Sr+k’h} 
> 


x(b) o(b) 


Siegel, 1, Lemma 27 and (81). 
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on putting 3 = x+y. Now 
givin: — b” when 26x = n (modb), 
(0) 0 otherwise. 
Hence 
G(S,r,b)|2 = b™ SY etritrr’ex+x'6) 
¥(b) 
2S¢=n(b) 
<b™ > 1=06"N(26,b) < b™2"8 
b) 
2e¢—nid) 
by (7). We must consider 
g2G,r) = > etminver (9) 
¥(b) 
25r=n(b) 
more closely. 
(i) Suppose 6 is odd. Then 2Gzx = n(modb) implies Gr =n 
(mod), so that, as rx’Sx is an integer for these values of xz, 
g(2S,r)= YS 1=N(G,d). 
x(b) 
Ex=n(b) 
(ii) Suppose b is even; say 6 = 28. Then 2Sz = n (mod 6) implies 
Sx =n (mod 8) so that, putting x = By+3, 
g(2S, ry a e2tire’ Sr 
x(2B) 
Sr=n(p) 
e27ir(By +3)’ S(By +3) 


3(B) (2) 
S3=n(B) 


e27irs'S3 eTiaBy’Sy 
if BS is odd, 
2™9(S,a/B) if BS is even, 


m 
{ 


> eTiaBy’ Sy _ i 


1 
> eTiaBsxeyz | 
0 j 


(2) c=1 'Ue= 


(0 if some Bs,,, is odd, 
| 2” if Bs,, is even for k = 1, 2,...,m. 


(ii) (2) Suppose f is odd and GS is even. Then (1/8)z’Sz is an even 
integer for all integers x such that Sx = n (modb). Hence 


9(S,a/B) = 
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(ii) (6) Suppose f is even. Putting x = By, as x runs over a com- 
plete remainder system to modulus f so does y. Also 


Sx = n (mod) 
is equivalent to Dy = n (mod), 
i.e. to €.Y, = O(modp) (k = 1,2,...,m). 


The solutions of this system of congruences which are incongruent 


to modulus f are given by 


y;, = (B/o,)z, (& = 1,2,...,m), 
where z, runs over a complete remainder system modo; With these 
transformations 
x’ Sx = y'B’SByn = yn’ WDy = f3'C3. 


Thus o1 Om 
g(S,a/B) = > ... > eree%s 
“4=1 lm=1 
m OK . 
> emicate), 
k=1‘x% 
rs m 
since a3'€3 = > Cp, 2, (mod 2). 
K=1 
Now, if c,; is odd, then », = 0 and p; = p, so that o, contains the 
same power of 2 as f, i.e. o;, is even. Thus 
a 0 when c¢,, is‘odd, 
SS eticnz — / kk 2 
ae \o, when c,; is even. 
Hence 
(0 when € is odd, 
| o,...0, = N(S,B) when € is even. 


m 


Q 


g(S, a/B) 
This completes the proof of the lemma. 
5. The function /(S, 7) 


We need the following inequalities for |f(r)|, where 


f(r) a f(S, rT) > emitr’ Sr 


x 


They will serve instead of the modular transformation equation in the 


Hardy-Littlewood treatment. 
Lemma 4. Suppose a, b are integers such that (a,b) = 1(b > 0). Let 


A= : 
~ {(2a—ba)?+b*y?} " 


w.= 2a—br, T= x+y, 
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Then f(S,7)| < |2/oo]m2f(S", ia) 
and, when G(S,a/b) $ 0, 


‘ie eae 1 /N(26,6)\3 : ae 
If(S,7)| > = wa = Mia y e-™™’s ‘| (y >0). (11) 
i | xr 
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(y > 0), (10) 


Proof. As all the series involved are absolutely convergent when 


0, we may rearrange the series by means of the substitution 
= by +3, and obtain 


x t) rm > evitx’Sr > > e7i{(2a—w)/b}(by +3) S(by +3) 
»y 3(d) 


aia > e27i(alb)3’S3 > e—Ti(w/bXby +3)’ S(by +3) 
3(b) v 


1-4 
i e27i(ald)s’Ss 8 aes Se il(I/bay’S—*y —(2/b)y’3} 
(thw)? : 

3(b) y 


by the reciprocal functional equation. That is 


AE) = Gay me 2, AE (S, a/b, (—1/b)y)erabow'e'y, (12) 


Consequently, by Lemma 3, 


f(S,7)| < |2/ea|™2 ¥ e-*9'S™"y — |2/eo|™2f(S-, 0) 
3) 
and so 


“- a |O(S, a/b)|—{b™N (2G, b)}! py adi: 
\bw |" 


i (N(2S,6)\3; 5 ibid’ 
anal S ai —_ 2 a ‘I, 


when G(S,a/b) # 0, since then by (5), (6), and Lemma 3 


|G(S,a/b)| = {b"N (26, b)}. 


6. Asymptotic behaviour of {(S,7) 


The behaviour of f(S,7) near the various classes of points on its 
natural boundary will now be discussed. First of all we have trivially 
TuHeorEM 1. Jf x is rational and equal to 2a/b where (a,b) = 1 
— iets 
» > 0, then 2 G(S, a/b) 
f(S,7) ees Stbmym? 
as y—> +0. 
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Proof. In (12) of Lemma 4 
w = 2a—b(x+iy) = —iby, 
so that 
f(S,7) oft > G(S, a/b, (—1/b)x)e lo*y)e’S~"e 
5,7) = lS, a/o, (— ila a © 
: (b2y)m2 - / / 
__ GS,a/b) 
pien Sipmyml? 1’ 


where by Lemma 3 


2 \m/2 ; as pa 
oat < 1 > e (7/b*y)x’'S “x — o(y~™?), 
by 
5 


We shall suppose from now on that 2 is irrational. Let 
l ] 


a,+ G+ 


ir = a)+ 


and p,/q, be the rth convergent to x. Write a;,, for the complete 


quotient corresponding to a,,, and 


, , ' 

Gr = Gi Wr Tv 
so that |4aq,—p,| = 1/9y44- 
We shall say that y lies in the range R, if 

9ly’2 < <— 9/2 

210r1 SY S 2/97 
i.e., putting y = 20/4, ria. 9 = Ur+alIr 

ijn <0 < 7. 
Since 9-41 —Grana = (4-41 —Gr4)9, > 9, 
when 2 is irrational, we see that the range R,,, overlaps R,. Hence 
as y tends to zero it will always lie in at least one range R,. 
Lemma 5. If x is irrational and y lies in the range R,, then 
f(S,7)| < K[min(q}?.,, 4/q2 y?)]}"4, 
where K is a constant depending only on S. 
Proof. Take a = p,, b = q, in Lemma 4 so that 
a a - 1/a’2, . -1¢2 42) 
jo! = |2p,—- Gr ®&— 19, Y| = 2CL/Grpa age y*)* 
4  -£ 
and A=—_, Y — : = ae 
4driatGry? = 1+e* 29, 

We notice that A has its least value in R, when 


f=, or I/y 
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Gr i 1 
2(q2-+ G1 


Hence f(S >~-1 , iA) “eo 4, 


so that in R, A> 


< Ye-twe's — K, 
® 


Thus by Lemma 4 
INS,7)| < K(/gF.+dary?) 


< K[min(q?2,,, 4/q2y?)}""" 
when y lies in R,. - 


THEOREM 2. For any irrational value of x 


f(S,a+iy) = o(y-™?) 
as y > 0. 


Proof. By Lemma 5 
If(S,a+ty)| < 2™2K(g,y)-™? = o(y-™*), 
since as y > 0, r > o and q, > ©. 
THEOREM 3. If 3x is an irrational number with bounded partial 


quotients, then f(S,x+iy) = O(y-™!4) 


as y > 
Proof. We have 
Grit = F414 G +1 < (G-42+2)q, = O4-) = 
when y lies in R,. Thus Lemma 5 gives 
IS, x+ty)| < Kg = O(y-™). 


TueoreM 4. Let (t) be any function which is positive, continuous, 
and monotonic increasing for large values of t and such that 


ie 


is convergent. Then for almost all* values of x 


renin GA) 


as y->0. In particular, for almost all values of x 
f(S,a+iy) = O(y-™*log?"*<(1/y)). 


* i.e. apart from a set of Lebesgue measure zero. 
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Proof. By Khintchine’s lemma,* for almost all values of x 


s0q—p| > (p,q) = 1, 


1 
q¢(q)’ 
for all but a finite number of values of p and g. Hence, 

1/Gr41 = |3@9,—P,| > 1/4, $(4) 
if r is sufficiently large, say r > ry. Thus, if y lies in R,, 
Gr+1/% < PG) < 4) (r > 1). 
Substituting in Lemma 5 we obtain 
ba ; 2 ,/1\ ys 
AS. 2-+iy)| < K] =a) ]" 
THEOREM 5. If x is an algebraic number of degree k 
f(S,a-+iy) = Oy-*) 
as y > 0, where x = m(k—1)/2k. 
Proof. As 4x is algebraic of degree k, we have by Liouville’s theorem 
1/¢¢4, = |$29q,—p,| > A/g*-}, 
where A is a positive number depending on x only. We divide R, into 
two parts. 
(i) If 2/¢,¢-41§ < y < 2/q?, we have 
Y > 2/4941 > 2A/q*, 
so that 1/q¢, = O(y"*). Hence 
If(S,w+1y)| << 2"?K(g,y)-™? = O(y~*). 
(ii) If 2/g/2., < y < 2/¢,9¢/4,, we have 
Grr = O(9, G11)" 4} = Oy"). 
Hence f(S,7) = O(gim'?) = O(y-*). 
By (i) and (ii) the result holds throughout the range R, and so 
generally. 
7. Q-theorems 
Let yw be the smallest number such that 
by e-men ls < 3 
F 
when ¢ > pv. Then yp exists and depends on S only. Thus in Lemma 
4, if G(S,a/b) A 0 andA > gp, 
If(S,7)| S FS ew|-m?, (14) 


* A. Khintchine, 1. 
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In particular, if g, is odd, we see from Lemma 3 that 


G(S, P,/qr) # 9. 
Now, if a = p,, b = q,, and y = 2/9,9/41, 
A = 941/49, jw| = 2V2/97 44. 
Hence we have 
Lemma 6. Jf y = 2/9,4/41> Ua1 > 4g, and q, is odd, then 
Lf(S, 7)| > 22-3" S-igimi2. (15) 
THEOREM 6. Let w(t) be any positive continuous function which 


tends monotonically to zero as too. Then corresponding to this 
function there are some irrational numbers x for which 


f(S,a+iy) = Q(y-™?aw(1/y)). 

Proof. It is enough to define a suitable sequence of partial 

quotients. First we notice that, since 
PraiG—-PGu = (— 1 y’, 

it follows that q,,, is odd when q, is even. Having defined a,, a, 
in an arbitrary fashion and hence p,/q,,..., 9;./q;.. we proceed to define 
the remaining partial quotients according to the following rule: 

(i) If q, is odd, choose a,,, so large that 


/2 


w(4a,.,9¢7)<q7™"" and a,,, > 4p. 


3 
(ii) If g, is even, choose a,.,, arbitrarily, e.g. a,,, = 1. 
In the latter case g,,, will be odd. Hence the number so defined is 


such that for an infinite sequence of values of r 


G41 > By +19p - 4uq,, 


oe ~ 2 te ; 
qo™t > w(4a,., 9) ae (39, %41) = of) 


and g, is odd. Hence by Lemma 6 
If(S, T) > 21—3m S-tg'm 2 


me 21—-5m/2,S—-by- mag m/2 
~ 21—5m/2 V-4 jo ] 
~ Zl—sm 2S—ty-mw ; 
y 


as y—>0 through the sequence of values corresponding to the 
sequence of values of r. The result follows. 
We next require a slightly modified form of Khintchine’s lemma. 
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Lemma 7. Suppose that x(t) is a positive, continuous, monotonically 
increasing function such that 


is 8) 


| dt 

typ(t) 

is divergent. Then for almost all values of 0 

] 
iga—v| < 
wg—P| < 
q¥(q) 
has an infinite number of integer solutions p, q such that 
(p,q) = 1, q>0 and q=1(mod2). 

The only difference between this and Khintchine’s lemma lies in 
the extra assertion that q is odd in an infinite number of the solu- 
tions. It is easy to adapt Khintchine’s later metrical proof* to 
justify this addition. 


lod 


THEOREM 7. Suppose that y(t) fulfils the conditions of Lemma 7 


and that b(t?) Ofb(t)} 


ast—> 0. Then for almost all values of x 


f(S,x-+iy) = a( *4(I »|") 
y 


In particular, for almost all values of x 
I(S,x+ty) = Qty" log”*(1/y)}. 
Proof. By Lemma 7 for almost all x 
\d0q—p| < 1/qh(q) 
has an infinity of solutions such that g is odd and (p,q) = 1. But, if 
|saq—p| < 1/29, 
then} p/q must be a convergent of $x. Hence for almost all x there 
exists an infinity of convergents with odd denominator such that 
V/Gp+1 =e \52q,—P, a 1/9, (q,). 
Let Q denote this sequence of values of g,. For all sufficiently large 
g, in Q we have 
=* H(g,) > Ay. 
Hence, if y = 2/q,q;,, and q, is sufficiently large in Q, Lemma 6 
atin f(S,7)| > Cami? > Cy-m*{p(q,)yn's. (16) 


* A. Khintchine, 2. + See, for instance, O. Perron, 1, 45. 
t The constant C is not the same number in each occurrence, but it is 
always independent of q, and y. 
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Now* for almost all x Ge41 = O(?). 
Hence for this set of numbers x 
b(1/y) < WO, Gre) < ¥(Ca?) < CH(q,), (17) 
provided that r is large enough. Combining (16) and (17) we obtain 
the desired result. 
8. A lattice-point problem 


Let R(N, 2x) om R(S, N,2x) _ a _e7izQy 


am -> A(S rjemizn, 


where A(S, 7) is the number of integer > y such that 

Oy) = n. 
The Q-theorems proved for f(S,7) easily lead to analogues for 
R(S, N,2z). 

THEOREM 8. Let w(t) be any positive continuous function which 
tends monotonically to zero as t- 0. Then there exist corresponding 
irrational numbers x for which 

R(S, N, x) = QIN™?w(N)}. 

Proof. Suppose the theorem false. Then for any irrational 

number x R(S, N,x) = o{N™%y(N)}. 


Thus for any irrational x 


f(S,7) = FE 4 A(S n)eman myn 


n 


= (l—e-7¥) ng R(n, x)e-7" 


rr 

=0 (y, > nme (n)e-7"), 

We may suppose without loss of aiding that 
w(t)>1/t as t>oo. 


Hence, putting Y = ih 
fE,7) = 0(y = + 0(yoo¥) 


n= n= 


x 
nmite-mvn) 
Pi 
) 


on o ( y m4 o(y- m2) (1/vy) 
on vd m2, oa (1/vy)) 
for any irrational x. But this contradicts Theorem 6, hence the 


result holds. 
* See Khintchine, 1. 
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THEOREM 9. Suppose that ib(t) is a positive, continuous, monotoni- 
cally increasing function such that 


x 


r at 
tip(t) 
is divergent and w(t?) = Ofeb(t)} 
ast— oo. Then for almost all x 
R(S, N, x) = Q([Np(N)}*) 
In particular, for almost all x 
R(S,N, x) = Q(N™*log™*N). 
Proof. Suppose the theorem false. Then for an 2-set of positive 
outer Lebesgue measure 
R(N,x) = o([Np(N)}™"*). 
Hence for an x-set of positive outer Lebesgue measure 
f(r) = (1—e-7") y R(n, x)e-7" 


_ 
n=0 


, x * 
= 0(y > [np(n)]™4e-7¥" ) 
‘ n=0 . 


: oly 


(t) a. 7, aryt a 


0 


= @O (y -ma [ | eeafs(au/y) J/4e-7™ an) 
0 
Now, as we may suppose #(t) > 1 fort > 0 
(uly) < maxf{y(w?), b(1/y?)} < b(u*)b(1/y?) 
= Ofys(u)p(1/y)} 
Hence for an x-set of positive outer Lebesgue measure 


f(r) = o(y-™4*{h(1/y)}™"*). 


But this contradicts Theorem 7, hence the desired result is true. 
We conclude by deducing 
THEOREM 10. Let d(t) be any positive, continuous, monotonic in- 
creasing function such that 3 


| td(t) 


as convergent. Then 


x 1  |m/2 
> {A(S, n)}?eP o(| 4(1/0)| ) 
n=0 y 
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In particular, for any « > 0, 


2) 


SAS, 2) Pe ‘as Ofy-m? log? +«(1/y)}. 


n= 


Proof. We have 


> {A(G, n)}2e-277" = 


n=0 


O(L(1/y)o(1/y)}"*) 

by Theorem 4. 
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AN APPLICATION OF FOURIER TRANSFORMS TO 
THE EVALUATION OF A CLASS OF REPEATED 
INTEGRALS ARISING IN CALCULATIONS OF 

PROBABILITY 
By J. A. CARROLL (Aberdeen) 
[Received 28 December 1937] 

1. A cuss of repeated integrals occurring in some physical problems 

is of the form 


| ‘a | fultw)fnaltn 1): fi(t,) dt, dt, ... dt, 


where the variables ¢,, ty,..., t, are restricted by the condition 





8 <t,+t,+...+4, < s+ds. 
This can be reduced immediately to the form 
ds | Sn(8 Yn fn Yn a—Yn—2) +++ fo(Yo—- YF) Uy ++ na, 


i.e. an n-tuple integral is reduced to an (n—1)tuple integral and 
the ranges determined. 

As an example of the origin of these forms, let f(x) dx be the 
probability of a displacement between x and x-+-dz, due, say, to 
Brownian movement in a certain interval, then, taking 


h=h=-=h=sh 
ds ae S(&_) f(%_~1) --- f(2) de, ...dz,, 
where s < >} a, < s+ds, is 4(s)ds, the probability of a displacement 
between s and s+ds occurring as the resultant of n successive dis- 
placements under like circumstances. 
Transforming, we have to evaluate 


[ eee ( f(8—Yn—DL£Yn—-1—-Yn-2) «+ Fa(Y2—-Y A(Yr) Wt «+ Yn a 


Analytically this expression is unwieldy and its properties, as a 
function of s, can only be extracted with difficulty. The numerical 
evaluation of ¢(s) when ” is more than 2 or 3 is clearly a very 
laborious task. 

It is the object of this note to show that this problem is capable 
of much simpler and more elegant analytical formulation, and 
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further that thereby the numerical computations referred to are 
greatly simplified even for a rigorous solution, and can be reduced 
to a very simple general approximate form applicable in nearly all 
practical cases. \ 

2. Let the functions f,(t) be bounded along the real axis, positive 
at all points of it, and tending to zero at infinity in both directions 
so that for each f, the Fourier transform g, exists and has the usual 
properties so that « 


= Fa0 je the-™ dt, 
1 q iu 
1) = Jo. J g(ue du. 


If f be not continuous over the whole range, the integrals are to exist 
as Lebesgue integrals. Further, let the functions be normalized so that 


— eo 


and the origin so chosen that 


@ 


[fo dt = 0. 


We now proceed to build up step by step the repeated integrals 
to be evaluated and show that the Fourier transform at each stage 
is simply the product of the transforms of the individual functions. 

Let F,(s) = fils), 

F(s)ds = | | fala) f(x) dey der, 


F,(s)ds = [ ... | fly) --- Fults) dy -..dity, 


where s < 4 pee .+a, < s+ds, for F,(s). Now consider Fj, where 
8 <4,+a, < s+ds, and put x, = t, x, = s—t, then 


8) = { fals—t) fx(t) at 


-2 


and it is readily seen that 


= f t6 —t) F_,(t) dt. 


3695.9 N 
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fe) 


| F(t)e-*4 dt 


— 


; v( 277) 


5) fe af k. (t—v) F_,(v) ded. 


By Fubini’s theorem we may change the order of integration and 
»btain ‘ . 

G,(u) = (2z)!g,(u) G,_4(u). 

Thus G,(u) = (2m)™Ng, (wt) gy, -4(te) ... gu (e) 


2 
and F,(s) = (27)i- [ ef" (U) Jn—y(U) ... 9,(u) du, 
—w 
where F\(s) is the required ‘distribution’ function. 
As it sometimes happens that the g’s are known (or easily found), 
this formula is itself a useful aid, and simplifies numerical work. 
However, we proceed to a still greater simplification. 


3. An asymptotic expression for F,(s). 

The method of steepest descents leads at once to the result that 
for large n the form of F,(s) is e-“* and determines a in terms of 
the second moments of the individual distribution functions /f. 

, 
We have bd (919e---Gn) = 9199---Gn S Ir 
du Ir 
and the least zero is at uw = 0, where there is in general a maximum 
maximorum of G(u). All g’s possess this same property and 
g(0) = 1/,/(27), in virtue of the normalizing conditions. 
Thus near u = 0 


9 


u~ 


: us ) 
oo os g Mat ogha fy 


where = He) t) dt. 


1 : —3p2 . 
Thus g(u) elogg(u) — * e—tveu rr 4a 5 M2 ust O(us)| 
(277) | 24 J 
(omitting the term in u® which disappears on subsequent integration). 
1 


rr 2 1 3 2) 4 ) 
s G —Mue'| yo a (Ha Hs | 5 
Thu »(u) any" (mn 7 +0O(u )» 
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n 
where M = > 2, and diverges to oo with n, and 
p=1 


i 2) 


pF2 — | tf, (t) dt. 


—@ 
@o 


1 


- : I . 
Thus F.(s) ~ se | e-tn-t dew — nO (IIT) 


2 4 a{(2ar 


i) 


This satisfies the condition of consistency, i.e. that { F(s) ds == I, 


and is the required result. 
In the commonly occurring special case where 


Fa = faa ies =f, =f say, 
M = mpy, 


1 - 
and F.(s)~ “5 a imps), 
V (272 


The order of the omitted terms will be 


(1 (4-343) 
(Sn w§ 
of the leading term, and this tends to zero as n tends to 00. 
Incidentally it may be remarked that this constitutes a proof that 
the resultant of a large number n of sources of errors of observation 
is distributed according to the ‘normal error law’, however the 
individual sources may be distributed. This is a well-known result ; 
this note adds to it by showing that the ‘convergence’ is rapid, and 
allowing the parameter of the distribution function to be evaluated 
and the order of approximation to be calculated, for relatively 


small n. 


4. An actuarial problem 

In practice these results are much more convenient than is at first 
apparent. For convenience of reference by those wishing to use the 
results without disentangling them from the proof, I cite the example 
of an actuarial problem. 

A typical question, of practical importance, requiring too elaborate 
computation for answering in the ordinary course is as follows. A set 
of individuals contract to pay an annual premium during life, from 
entry, to a pension scheme until a fixed age, say 65, and on survival 
to that age enjoy a simple life annuity of appropriate amount. The 
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ordinary computations of the ‘net’ premiums for an annuity (that 
is, excluding the ‘loading’ of the premium by adding an amount 
deemed sufficient by the insurers to cover their operating costs and 
provide a suitable profit) are based on life tables and are such that 
the expectation of either a profit or a loss is just one-half. It is 
clearly important to assess the probability that the profit or loss 
on a group of » individuals will exceed some specified figure. This 
problem is essentially of the type dealt with above. 

For simplicity of illustration we will consider a group of n in- 
dividuals, all entering simultaneously at age 35 on a scheme to pro- 
vide at 65 a deferred annuity of amount A per annum for a premium 
(not returnable) of » during the 30 years’ interval of deferment. If 
5, represents the probability of dying at age x, and v the present 
value of £1 due in one year at the rate of interest used in calculating 
p {= AN,,,,,/(N,—N,+,)}, then considering a single individual we 
have, reckoning in terms of present values at the inception of the 


scheme, 
5, is the probability of a profit p 
Sn41 ~ . »  pi+v) 
]—y"t1 
xrtn ” ” = p(itve+...+v") =p ; 
—Wv 
] ‘a yn +1 
O4041 2” Pr) 9 p a au prt 
1—y"+1 41-07 
Srintr 29 ” ” Pp pe — Ayntt —." 


We have our usual conditions 
+3, =1, > (profit)s = 0. 
We require 1, the second moment 
= > (profit). 


There will as a rule be two ways of gaining a given profit P (either 
by death before maturity or death soon after maturity), but it is not 
necessary to find the resultant of these as in > (profit)*5 the relevant 
term will be P?(5,+6,), i.e. we need only form the second moment 
HM, of the table as it stands, a simple computation. The required 
information is then given by 

1 


P(e) = et, 
(272%) 
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where P(s)ds is the probability of a profit (or loss) between s and 
s+ds, and the likelihood of a loss exceeding L is 


o 


1 | e—(s"/2H2n) ds, 


4 (27p42%) 
L 


and reference to a table of error integrals gives the result immediately. 
To quote a numerical example, based on life tables and commuta- 
tion columns at 3} per cent., we find: 


For an annuity of £800 per annum the net premium is £95 per 


annum, and [lg = 32677, pg = 21-6144, 


taking of course A = 1 in normalizing the distribution function. 
Thus the probability of a profit or loss s to s+-ds is 


1 


e-(8'/6'5354n) dg. 
4 (6°53547n) 


the neglected term being 1/8 of the leading term; thus for n = 30 
the error is 0-4 per cent., for n = 3,000 it is 0-004 per cent. The 
table below shows the use of the expression obtained. 

If all the members are not of the same age and the time chosen 
is not at the inception of the scheme, the calculations are only 
modified in the relatively trivial respect that », must be found for 
each member, and M, the sum of these, used instead of n times po. 


Balance Sheet for Deferred Annuity Scheme at inception, with n 
members all at age 35 contracting for an annuity of £800 p.a., at 65 


Mean present | Present value of loss 
value of future | sustainable will exceed Additional pre- 
Number premiums or | figures below once in mium required to 
of annuities 4 times 10 times Net cover risk of loss 
members payable (a) (d) premium (a) (b) 


£ £ £ 3 £ £ 
30 47,530 5,344 10,153 9% 10-7 20-3 
120 190,120 10,688 20,306 9! 5:3 10-2 
3,000 4,753,000 53,440 101,530 9! . 2-0 
12,000 19,012,000 106,880 203,010 95 0-5 0-1 





The risk contemplated here is that of an unusual proportion of the annuitants 
surviving to a great age and so throwing a large burden on the scheme. The actual 
sum lost can be computed, but for simplicity it is here stated as its present value 
at the inception of the scheme. The corresponding value at the start of payment 
of annuities is approximately 2-8 times these figures and for the time (60 years later) 
when all members can be supposed dead 2-8 times this again. The importance of 
a large number of participants is clearly seen. 














THE MINIMUM MODULUS OF INTEGRAL 
FUNCTIONS OF FINITE ORDER 


A. J. MACINTYRE (Aberdeen) 
(Received 12 April 1938) 
In Pélya’s theory of gap theorems an important part is played by the 
following 

THEOREM. If f(z) is an integral function satisfying log M(r) = o(r), 
where M(r) is the maximum modulus of f(z) for |\z| =r, and «a > 0, 
then |f(n)| > e-% for almost all integers n.* 

As many difficult proofs of this theorem and its generalizations 
have been published, it may be worth showing that these results are 
trivial corollaries of Hadamard’s minimum-modulus theorem, at least 
when the precise form to be derived from a well-known theorem 
of Cartan is used. Cartan’s theorem states: if H is any positive 
number, the polynomial P(z) = (z—a,)(z—a,)...(z—a,,) satisfies 
|P(z)| > (H/e)" for z outside a system of at most n circles the sum of 
whose radii is at most 2H. We also require Nevanlinna’s formulat 





lf | Rr? 
[f(z)| =— | log|f(Re®)|- 
log} f(z)| an | og | f( Re )| Ra BFR coslO—) dO + 
R(z—a,,) 


+ > log) ———_— (R>r; z= re'?) 
R?—2a 
lan|<R " 

* The theorem was stated by Faber, Jahresb. der Deutschen Math. Ver. 16 
(1907), 285-98. A proof was indicated by Pélya, Comptes rendus, 184 (1927), 
579-81. For generalizations and further references see Pennycuick, J.of London 
Math. Soc. 12 (1937), 267-72. The applications to gap theorems are developed 
by Pélya, Annals of Math. (2) 34 (1933), 731-77, in particular 745-54. 

+ H. Cartan, Ann. Sc. de? Ecole Norm. (3) 45 (1928), 255-346. The theorem 
quoted is proved in pp. 272-6. There is an omission in the proof which is 
rectified in the version given by Valiron, Directions de Borel des fonctions méro- 
morphes, Paris (1938), 11-12. For Theorem 2 we only need a much earlier 
result due to Boutroux, Acta Math. 28 (1904), 97-224 (129-33) and Theorem 2 
itself is practically equivalent to Theorem 27, p. 89 of Valiron’s Lectures on 
the general theory of integral functions (1923). A. Bloch, Ann. Sc. de l Ecole 
Norm. (3) 43 (1926), 309-62, states without proof a theorem very similar to 
that of Cartan. Both Bloch and Cartan proceeded to other matters without 
explicitly stating the inequality (1), which is in effect the conjunction of 
Bloch’s lemmas 4 and 5 (pp. 320-3). Formal theorems on the minimum 
modulus of a function regular in a circle arising out of Cartan’s theorem have 
been given by H. Milloux, Proc. Phys. Math. Soc. Japan (3) 12 (1930), and 
V. Bernstein, Annali di Mat. (4) 12 (1934), 179-81. (For the last two refer- 
ences I have to thank a referee.) 

t R. Nevanlinna, Le théoréme de Picard-Borel et la théorie des fonctions méro- 
morphes, Paris (1929), p. 24. 
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where f(z) is supposed an integral function and a, are its zeros. 
Evidently 

R—r 
log| f(z)| > — 


R tbs 
m(R,f) — oeee (R, Uf) + > log| et, 
|a,|<R 


2a 
where (r,f) = = | log| f(re®)| dé 
7 
0 


From Cartan’s theorem with H = hR we obtain 


log] f(z)| >a m(R,f) — FA | R,Uf)—n(R)log = (1) 
outside a set of at most »(#) circles the sum of whose radii is at most 
2hR, where n(R) denotes the number of zeros of f(z) of modulus less 
than R). If h is small, it is evident that this inequality is satisfied 
for the majority of values of z, but there is some interest in deter- 
mining the largest regions of various types in which the inequality 
holds throughout. The natural regions to consider are circles and 
annuli of the form a < |z| < 6. Let us draw corresponding to each 
circle of radius h, R in which (1) is not true a concentric circle of 
radius h, R+86Rn(R)-*. The area included in this family of circles 


is at most 

aR? > {h,,+-5n(R)+*} < wR*{4h?+ 4hdn( R)-*+-8?}. 
Points outside these circles are the centres of circles of radii 8Rn(R)-* 
throughout which (inside |z| = R) (1) holds. If h and 6 are small, 
this applies to the majority of values of z. In order to obtain a more 
simple result let us suppose that f(z) is of minimum type with respect 
to a proximate order p(r), that is 


log M(r) = of{V(r)}, V(r) = re), 
lim te =k (0<k<o), 
Tmo r 


and let us confine our attention to values of r between }R and $R. 
Then, if h, 5, and « are positive, no matter how small, (1) gives 


log| f(z)| > —aV(r), 


provided that r (and hence R) is sufficiently large. Also the radius 
5Rn(R)-+ > rV(r)- for sufficiently large 7, and we have 
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THEOREM 1. Jf f(z) is of minimum type with respect to a proximate 
order p(r), then for any positive « 

log| f(z)| > —aV(r) 

throughout the circle |jz—l| < rV(r)-* for almost all ¢. 

Similarly, if we associate with a circle in which (1) does not hold, of 
radius h, and centre distant r,, from the origin, the annulus 

r,—h, R—SRn(R) < |z| < 7,+h, R+6Rn(R)-, 


the total measure of such values of |z| is not greater than (4h-+ 28) R. 


Arguing as before, we have 
THEOREM 2. If f(z) is of minimum type with respect to a proximate 
order p(r), then for every positive « 


log| f(z)| > —aV(r) 


throughout the annulus r—r/V(r) < |z| < r+r/V(r) for almost all r.* 


/ 


It is evident that this result includes the result of Pélya by taking 
p(r) = 1. 


* In the special case V(r) = r? these two theorems become the Theorems 
4 and 5 given by Pennycuick, loc. cit. 

The method extends to functions regular in an angle by using the corre- 
sponding form of Nevanlinna’s formula. See Nevanlinna, loc. cit. 1-3 and 
Acta Soc. Sc. Fennicae, 50, No. 12 (1925), 29. It may be that referred to by 
Pélya, Annals of Math. loc. cit. 745 (in the footnote). 














ON THE OCCURRENCE OF MILNE’S SYSTEMS OF 
PARTICLES IN GENERAL RELATIVITY 


C. GILBERT (Dundee) 
[Received 13 February 1938] 


1. THE systems of particles constructed by Milne* are of two types. 
These are (a) the simple kinematic system comprising a set of funda- 
mental particles having a hydrodynamic motion, (b) the statistical 
systems defined by certain space-velocity distribution functions. A 
universe consists of the system (a) and a system (b). The system (a) 
is then a representation of a universe in which the mass of all the 
particles of the statistical component vanishes, and has therefore 
been called the substratum or smoothed-out universe. The types of 
equivalence which these particles possess have been clearly stated 
by Milne,t and I only wish to recall that both systems are described 
by a set of observers moving with the fundamental particles of the 
system (a), and that in the experience of every one of these observers 
the same picture of world events is obtained. 

The first step in attempting to describe these systems according 
to the theory of general relativity is to determine the metric of 
space-time for the universe, and, since time and distance measure- 
ments have a significance only when they are taken in relation to 
the fundamental observers, the metric must be one corresponding to 
the system (a). The most general form of metric which, according 
to general relativity, corresponds to this system, has been shown by 
Walkert to be 

ds? = (R(x) {ae— alda?+dy?+-de%). (1) 


In this space-time the light tracks are given by the null geodesics, 
and the time and distance measurements are made according to 
Milne’s light-signalling conventions from observations of clocks 
carried by the fundamental observers. The form of F(X) can be 
further determined from the condition that the geodesic equations 
shall correspond to the equations of motion of certain of Milne’s 
particles. 

* Cf. Relativity, Gravitation, and World Structure (1935). Subsequently 
referred to as W.S. When other references to papers by E. A. Milne or 


A. G. Walker recur, they will be signified by name and date. 
+ Proc. Royal Soc. A, 156 (1936), 62. t M.N. 95 (1935), 263. 
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2. According to Milne’s theory the acceleration of a particle of 
the universe, having the coordinates and velocity P, ¢, V, in the 
experience of a fundamental observer, is* 

1d/V\_ G() Pp VZ 
Yidt\Y:}) xX oy 
where 
P y2 Z 
XxX =?——, Y =l1——, : —> = —., 
o A : é XY 

The particle may be a free test-particle projected into space in 
any manner from one of the fundamental particles, and then Milne 
has showny that, in the presence of the substratum alone, the 


(3) 


acceleration function G(é) has the value 
@(é) = —1. (4) 
On the other hand, if the particle be a member of a statistical 
system having the space-velocity distribution 
b(E) 
c&Xiy! 
the acceleration function takes the form 
G(é) = —1— “S (6) 
(€—1)'h(€) 


where C is a constant, which has been shown{ to be a measure of 


dxdydzdudvdw, (5) 


the gravitational mass of the particles. 

3. In a space-time having the metric (1) the geodesic equations 
can be written in the form (2), if instead of G(é) we write G(X), where 
2X F’(X) 

F(X) 


There is therefore a correspondence between the trajectories of 


G(X) = (7) 


particles in Milne’s theory and the geodesic paths of general relativity 
when the acceleration function has a constant value.§ For example, 
the trajectories having G(é) = —2n correspond to the geodesic 
paths in a space-time having the metric (1) with 


“ {2 n 
F(X) -(3) ; (8) 
The constant ¢), having the dimensions of time, has been introduced 


ilne (1936). + Quart. J. of Math. (Oxford), 8 (1937), 22. 


M 
Milne (1936). § Walker (1935). 


na 
+ 
+ 
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into (8) to give s the dimensions of time. If we take n = 4, the 
metric has the form 


2 
dr? = (de — 5 (da? + dy*+-az")), (9) 
and the acceleration function has the value (4) for a free particle in 
the substratum. We shall subsequently show that this is the metric 
of space-time in which a description of Milne’s systems of particles 
can be made. We therefore call (9) the metric of the substratum. 

The following work is concerned with the construction of statistical 
systems of particles, when the observer’s space-time has these pos- 
sible forms, and is based on a recent paper by Walker* in which he 
derives the Boltzmann equations in general relativity. 

4. The metric given by (1) and (8) is conformal to the metric of 
special relativity defined by the fundamental tensor g;; having Gali- 
lean values. The fundamental tensors of the two spaces will be in 


the relation X 
ae a Pe los(3 (10) 
0 


For convenience we shall often use a suffix notation, the suffixes 
taking the values 0 and 1, 2, 3 for the time-like and space-like com- 


ponents of a vector. Then 
dx! dx* 


X = gy, 2"x*, Y= Ink ede’ 


Z= mo (11) 
and the metric (1) becomes 

ds* = g,;dx'da’. (12) 
Hence, if A’ is the unit vector tangent to the path of a particle at an 
event P, we have Gi, =1. (13) 
Then the element of three-dimensional solid angle dQ, determined 
by the range of directions (A‘,A‘+dA‘) relative to an observer at P, 


is given byt dQ. = (—G)*(Ay) dA dA2d8. (14) 


da! dx? dx 
Opp AE lt, B 1 11), and 
Wt? ad’ a Then from (10), (11), an 


(12), the components (A°, A1, A?, A®) can be written 


X\" 1 X\" u X\" v X\" w 
aj) Y? \@) yr @) Y? ) Yt) 


* Proc. Edinburgh Math. Soc. 4 (1936), 238. + Walker (1936). 


We shall write (u,v, w) for ( 
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Hence 
p @\* 1 a(A}, A2, A%) X\3" dudvdw 
- ee |. | ea Seek LE hs = 
atleast (z) Y? (u,v, w) (7) yi 
Therefore dQ =. nian =r (15) 


{1—(u?+-v?-+ w?) /c?} 
This is an invariant expression for all observers, and a constant 
multiple of dQ (actually B/c*) gives the velocity distribution of 
fundamental particles occurring in Milne’s theory.* 
The element of volume of a three-dimensional space orthogonal to 


is 2 ta we r 
d' at P, is given byt 


2\3n ° » 
dV = (—g)!A° daldx*dz? = (%) — (16) 


Hence a system of fundamental particles having the number density 
v and velocity distribution V = P/t¢ will have a space distribution 


a 2 3n ft 1. * iz eS 
vaV = (2) — (7) 

This gives Milne’s distribution of fundamental particles 
as miinigts — (18) 

c8{t?— (a? +-y? +2?) /c?}* 
: B é 3(n—}) 
provided that v= a ( | : (19) 
cts \ to 


In the particular case n = 3, the particle density has the value 
B/c®t3. In the space-time of the substratum, therefore, the density- 
distribution of fundamental particles is stationary and homogeneous. { 
This result has been obtained by Milne.§ 

5. A statistical system of particles is defined by a velocity- 
distribution function ¢(x, A) such that the number of particles cross- 
ing the volume element dV whose paths issue in the angle dQ is 


given by|| dN = d(x, d)dVdQ. (20) 


In order that particle number shall be conserved, 4(x,A) must satisfy 
the generalized Boltzmann equation 


i 

— At —_—— I, NIA* = 0,7 21 

Cart ant # (21) 
* W.S. 87. + Walker (1936). t W.S. 64. 


§ Proc. Royal Soc. A, 159 (1937), 171. i| Walker (1936). 
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where the Christoffel panic are “es by* 


=; do 
ie = = li, +8) a i= 5 In 5 (22) 
and the symbols Ij, for a flat space are zero. 
We define dimensionless quantities x, 4 by 
hk hk 
= 4 = ——, a= i dX _ Inz vr" (23) 
t t 2t, ds ty 


and we look for solutions of (21) of the form 
$ = c-% *(y, 6). (24) 
This function ® must be dimensionless in order to make the right- 


hand side of (20) a pure number. 
From (10), (22), (23) we find 


Pay ; 9 2n-1 | 
Ti, var = \ae ae xi, (25) 
~ Xto 
® yi om 1 26 0 +X 2n ) 
Ox" . 


2 iS 


i 8 a roe ab 
r ANAK = —| ——— + 2ny*")|— 
ani -( x To 5 


Hence equation (21) reduces to 


207 + ja a (26) 
| x 06 
This has solutions of the form 
6 1 # 
® = f(T) =I say, (27) 
xX xX 
where, with the notation of § 2 
6 te 2n-1 
aan & t= (3) ; (28) 


Then, from (15), (16), (24), the number distribution (20) becomes 


. dxdydzdudvdw 
NV = : ; 29 
aN = fl) aacx aT (29) 


We can now find the statistical distributions which are independent 
of the arbitrary unit ¢). Since f(¢) is dimensionless, (29) will, in 


* Cf. L. P. Eisenhart, Riemannian Geometry (1926), 89. 
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general, only be independent of ¢, when f(f) = C,,¢-!, where C,, is 
a dimensionless constant. We then have 


iN — C,  dadydzdudvdw 


c(é—1)' xiv? ins 


Comparing this with (5) we see that, in general, the statistical sys- 
tems of general relativity, in space-time having a metric given by 
(1) and (8), correspond to the systems of Milne’s theory having 


constant 
¢—1F° 

In the particular case n = $, however, the number distribution 
(29) can take the general form (5) of the distribution in Milne’s theory. 
This gives the distribution of free particles in the space-time of the 
substratum. 

We therefore conclude that, in a space-time in which a set of 
fundamental particles having a hydrodynamic motion can be con- 
structed, it is also possible to construct statistical systems of particles 
defined by the distribution functions (31), except in the particular 
case when the metric of space-time has the form (9). In the latter 
case the form of 4(€) is not uniquely determined. We shall now show 
that constrained systems of particles, having the properties of Milne’s 
statistical systems, can be constructed in the space-time of the sub- 
stratum. The argument is essentially similar to that given by Milne.* 

6. The particles of the statistical systems just considered were 
following geodesic paths, and in the space-time of the substratum it 
can be shown that € = constant along a geodesic path. We now 
suppose that each particle of a statistical system, in this particular 
space-time, is constrained to follow some path along which € varies. 
We therefore introduce into the geodesic equations an acceleration 


¥(€) = (31) 


four-vector F’, such that 
dx 


dr 


+ Pi, = Fi, (32) 


where Ti, \/A* is given by (25) when n = } and 2! = dz‘ /dr. From 
(9), (23), (25) this equation reduces to 


x! d{xi =, ( i 24 a awe 


Y at/e 


Yidt\¥? dt 


* W.S. 179. 
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This can also be written 


1 d/1 dz* , ZSdx\l_ GF 
nalrs a) + (2— Y ax i (33) 

The form of F*‘ must now be chosen so that (33) satisfies the 
cosmological principle, that is to say, is invariant under a trans- 
formation from the coordinate system of one observer to that of any 
other observer. According to Milne, equation (2) gives the general 
form of the equation of motion of a particle satisfying this principle. 
This equation can be written in the form 


1 aya) +(P— —_ p+4@(P— > ug (34) 








Y! dt\Y? Y /X Y/X 
We also associate with (34) the time-component equation 
1 d[c cZ cZ\ 1 
— — | — z +4 —=)|s- 5 
rian) t(¢-F)x -U+e@-F)y- 


Comparing these equations with (33), we see that 


pL QE)(_2 ae 


rR ea 
Using (9) and (11) this can be written 


Fi- [1 + G(é)] (fe NAMx!i 7 Fim xe") 
te x - 





= [1+ @(g)]2n% amai—ainm) (36) 


7. We can now construct the constrained systems of particles 
having the space-velocity distributions (5), in the space-time of the 
substratum. The distribution function %(€) must satisfy a certain 
generalized Boltzmann equation in order that particle number shall 
be conserved. When the equation of motion of the particles is given 


by (32), this equation has been found to be* 


OME) ye _ SA) sae + SHEP} = 0. (37) 
From (32) this can be written 
dyp(E) F‘ 
, 8 
dr +48) x ris (38) 


* This equation was obtained for me by Dr. T. G. Cowling. It can be derived 
by the method used in Walker’s paper (1936), on the Boltzmann equation, if 
account be taken of the additional term F"' in the equations of motion. 
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For the space-time of the substratum 
So Tea 0& — 2945 7"Inx atl. 
x* ane x* 
so that from (36) 


oF 291m Nam dG(€) i as . . 
oF’ — —MmMeraG16)¢_ yy 401+ 0063] 
Also, since (33) reduces to (2), it can be shown that 
dé 22,, 
= —x & 1)[1+ G(€)], 
Pern UE) ——_— 2D NX™ APE) |p ’ 
whence Fai dé (€—1)[1+ G(€)]. 


From (39) and (40) equation (38) becomes 





29m” 
x 
x | BP EDO 6O}+Heae—DTE + wore] =o 
dé dg 
1 dp(é) 4 age) ee (41) 
w(f) d&  14+G(€) d& © 2(g—1) 
This equation integrates to give Milne’s relation (6). 

8. We conclude that the substratum can be completely described 
by the space-time (9), in accordance with general relativity. It is 
also worth pointing out that our work has shown that a system of 
free particles having the general distribution function #(£) can only 
be constructed in a space-time having the metric (9). This is in 
agreement with Milne’s conclusion that G(é) = —1 defines a pure sub- 
stratum. We also observe that Walker’s result concerning particle- 
distribution in general relativity spaces agrees, as a special case, with 
Milne’s general relation between G(é) and #(€) 

The metric (9), besides defining the space-time of a pure sub- 
stratum, gives the interval dr for a universe. In the latter case, 
however, the tracks of particles are not given by the variational 
principle 5 { dr = 0, but depend on a constant C which is a measure 


Therefore 


of the mass of the particles. 
The interval of proper time of an observer in the universe is, 


from (9), dr = +t, dt/t. (42) 








ON MILNE’S SYSTEMS AND GENERAL RELATIVITY 193 


Here ¢ is the ‘kinematic time’ of the observer which we shall assume 
to be increasing with 7, and to have the value ¢, when +t = fy. 
Then, integrating (42), we find 

t—t, = th logt/ty. (43) 
This is the relation Milne has found* between the ‘kinematic time’ 
and the ‘Newtonian time’ of an observer, when ¢, is the time at 
which the clocks graduated on these different scales agree. The latter 
has been called the ‘present time’. The ‘Newtonian time’ of an 
observer is therefore the same as his ‘proper time’, and the constant 
i) is a measure of the present age of the universe on the kinematic 
scale. 

9. Summary. A description is given of Milne’s systems of 
particles according to general relativity. The metric of space-time is 
one for which the geodesic paths give the trajectories of free particles 
in the substratum. The acceleration function for the particles of a 
statistical system has the same form as that given by the kinematic 


treatment. The “proper time’ of an observer agrees with his “‘New- 


tonian time’. 
* Proc. Royal Soc. A, 158 (1937), 324. 

















ON DIRICHLET SERIES WHICH SATISFY 
A CERTAIN FUNCTIONAL EQUATION 
By H. HEILBRONN (Cambridge) 

[Received 12 January 1938] 


WE consider a Dirichlet series 
d(s) = Da,n- (1) 
n=1 


which satisfies the following four conditions: 
(i) the series (1) has a half-plane of convergence, so that 
x 
, . 
fly) 74,.¢"" 
n=1 
is absolutely convergent for y > 0, » > 0; 
(ii) 4(s) is regular for all finite values of s apart from a finite 
number of poles; 


(iii) we can find a positive integer k and positive numbers A and 


A such that (8) = A-*T*(s/A)d(s) 
satisfies O(s) = O(1—s); (2) 
(iv) $(s) = O(|t|”) (3) 


uniformly for o > 3, |t| >00, where the constant P depends on the 
function only.* 

By ®,(s) we denote the rational function which vanishes at infinity 
and makes ®(s)—®,(s) an integral function; 2,..., x, are positive 
variables; we denote their sum by S(x) and their product by N(x); 
p is an abbreviation for A/k. 

Then we have for all values of s 

r dx, dx 
, , §N (9)slA T(y)-siAr f Var) 1 k 
O(s) = ®,(s)-4 [ Ses | {N (x)** + N (x)*-"} f {AX S(x)} —... —. 
2 Ly 
“ ° 1 k 
N(x) 1 (4) 


The proof of (4) is easy. It follows from (3) that the difference of 
the two sides of (4) is regular for all s and bounded if o is bounded. 
Hence it suffices to show that the difference tends to zero as o tends 
to +co and —oo. But (2) shows that all the three terms in (4) are 

* This condition can be replaced by a weaker one by the use of an argument 


of Phragmén-Lindelof. 
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even functions of s—}. Hence we need investigate the case a > +00 


only. 
For sufficiently large positive o the series in (1) is absolutely con- 


vergent. Hence we have for n > 0 


A-*kT(s A)n-slk a | xslre-- Apt OH da 
x’ 
0 
A-T*(s/A)n— = e& [ N(x x) Ne—- nu g(a) OX) oe 
a J Z, Ly 
N(x)>0 
4 dz, da, 
> 4 s/A LQ» k 
(s) = & | N (x) "f,fAMS()} Tn 


Nix)>o 
sut the two integrals 


= [%e Af {AM S(x “a 


; z, Ly 
0<N(x)<1 
- dz, dz, 
and [ as | N (x)t-9AF {A¥S(2x)} —... —t 
J ; Z, Lp 
N(x)>1 


tend to zero as a > 00. So does ®,(s), and (4) is proved. 
We mention one application of (4). 
Let d,, d,, ds be discriminants of quadratic fields and let L,(s) 


denote the series : 
> (d;/n)n 
n=1 
Suppose further that d,d,d, is a perfect square. Then 


3 
) TI La(s) 
i=1 
has non-negative coefficients and satisfies our conditions (i)—(iv) with 
3 . 
k = 4(5+ > signd,), 
i=1 
A= hk, 
kn2 7 
A = 24-*y* TT |d,|-*; 
i=1 
and (4) is the identity from which Siegel* derived his important 
‘esu 
result log L,(1) = o(log|d)). 


* Acta Arithmetica, 1, 83—6. 

















ON CERTAIN HANKEL TRANSFORMS 
By A. ERDELYI (Brno) 
[Received 28 January 1938] 


1. Inarecent paper* R. 8. Varma has developed the Hankel trans- 
form of the function e~**x!’ L(x) in a series of Laguerre polynomials 


LY (x) = ~ ( 2 (—2x)* 


n—m] m! 


m=0 

and proved Wilson’s integral equation} to be a special case of his 
formulae. It is the purpose of the present note to point out that 
term-by-term computation of the Hankel transform of any function 
of the form e-**x# P(x), P(x) being any polynomial, yields a finite 
series of Whittaker’s confluent hypergeometric functions. If 2u = v 
(the order of the Hankel transform), this expansion reduces to a 
series of Laguerre polynomials. In the special case 

P(z)= G(r) and Yo=v=—=-« 
the series obtained by this can be summed by means of the multi- 
plication formula of Laguerre polynomials. The integral equations 
satisfied by Laguerre polynomials and products of such polynomials 


are special cases of this formula. 
. n 
2. We denote by P(z) = > o,,.” (1) 


any polynomial, and compute the Hankel transform of order v of 
the function . " 
F(x) = ea" P(x), (2) 


i.e. the infinite integral 
T = | J,{2J(xy)} Fy) dy, (3) 
0 


putting (2) in (3) and integrating term by term. Thus we obtain 
J,{2,/(xy)}e*uye+™ dy. (4) 


* R.S. Varma, Acta pontif. Acad. Sci. 1 (1937), 37-41. 
{+ B. M. Wilson, Messenger of Math. 53 (1923/4), 157-60. 
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Using Hankel’s generalization of “Weber’s first exponential integral’* 


in the form 


oe PA+P+)) 4 4-wanae (2 
| ekuyrJ,{2,/(ay)} dy = T'(v+1) saath Ze ”Myas(z) 


[R(k) > 0, RA+4y) > —1], 


we obtain at once the expansion 


0 


[ ety Pd {2y ley)} dy ) 
0 
1 ~~ ad 
a ke-# hate Halk) z Cm U(u+4dv+m+1)k~ ama] 


Pv+1) m=0 


[R(k) > 0, R(u+4v) > —1]. 


This expansion is also valid for power series P(x) convergent for all 
finite values of |x|, supposing that the integral on the left and the 
infinite series on the right converge absolutely. 

3. In (5) we put » = 4v and observe that in this case Whittaker’s 
functions can be expressed in terms of Laguerre polynomials, for 


m! T(v+-1) 
I'(m+v+ 1) 


zw Hee L)(z), 


M,,, m +3,4v(2) = 


Thus we obtain 


. n ® 
[etry PY), (2y ey) dy = k-r-Aahe-2 > meg k-mLS0(7). 


[R(k) > 0, Riv) > —1]. (6) 


The expansion of the Hankel transform in a series of Whittaker 
functions and Laguerre polynomials respectively is intelligible by 
the fact that these functions are characteristic functions of the 
Hankel kernel, the corresponding characteristic numbers (the only 
ones of the Hankel kernel) being +1, each of infinite order. 

4. Some special cases, in which the series on the right of (5) and 
(6) can be summed by means of known relations, may be considered. 
First put 

P(x) = Le (ke), ey = Mm PHT EMT pom, 
m!n! T(u+4v+m-+1) 

* G. N. Watson, Bessel Functions (1922), § 13.3 (3). See also A. Erdélyi, 

Math. Annalen, 113 (1936), 357-62. 
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and make use of the relation* 


* ~ 
— (—2)m 3 | a 
> m! My, +m+av() ~ T(vt+n- i M+4n+4,4v+in(2)> 
m=0 | 1 
thus obtaining 
o 


+ e-kuytJ,{2, (ay) Let (key) dy 


a ie ly ” 
_ P(u+ VENT) py intgdn—te—Hak) VU i % 

—«dP@w+n+)) ~ 

[R(k) > 0, R(u+4r) > —1], 
an integral representation of Whittaker’s function, which is equi- 
valent to a generating function of the confluent hypergeometric 
function.t Next replace x and y by kx and y/k respectively and 
put in (6) 
n+a\ (—)™ 
Pte) = Lee), oy = (PTC, yaa 


n—m} m! 


Using the multiplication formula of Laguerre polynomials{ 


> (252) —oen@@) = Ene *), 
, \v—m 1—k 


m 


we get at once the relation 


Ew z 
p-UzpaJ [2 ( x) ] riXe (a) > 
i eYyi~J {2,/ (ay) } LK ({) ay ( —;)" e*L' (; i (8) 
[R(a) > —1]. 


Wilson’s integral equation is the iailal case k = 4 of this formula. 
The limiting cases k > 0 and k > 1 of (8) yielding Le Roy’s integral 
representation of Laguerre polynomials$ and its inversion, the 
formula for the Hankel transform of the function e-*a!”L(x) re- 
spectively, are also known. 

* A. Erdélyi, Math. Zeits. 42 (1937), 641-70 equation (1, 1). 


+ A. Erdélyi, Monatsh. f. Math. u. Phys. 46 (1937), 132-56 equation (7, 3). 
t A. Erdélyi, Math. Zeits. 42 (1936), 125-43 equation (5, 4). 


§ E. Le Roy, Annales de Toulouse (2) 2 (1900), 317-430. Cf. especially 
379-84. See also G. Szegé, Math. Zeits. 25 (1926), 87-115. 





ON WARING’S PROBLEM 


By LOO-KENG HUA (Cambridge) 
[Received 1 February 1938] 


THE object of the present paper is to give a proof that Hardy and 
Littlewood’s asymptotic formula for the number of solutions of the 
Diophantine equation 
N = ak+...+a% (2, > 0) 

is true for s > 2*+-1. This result is new only for k < 14. Fork > 14, 
Vinogradow’s contribution is much better than this. The most in- 
teresting particular case is k = 4, for which Estermanny and 
Davenport and Heilbronn{t proved that every sufficiently large 
integer is a sum of 17 fourth powers, but they gave no asymptotic 
formula for the number of solutions. 

More precisely, what I am going to prove is the following more 
general theorem: 


Let P,(2) P(x) be s integral-valued polynomials of the kth degree 
and let their first coefficients be positive numbers a,,..., a, respectively. 
Let r(N) be the number of solutions of the Diophantine equation 

N = PB(2,)+...4+-B(z,) (x, > 0). 


Then, for s > 2*+-1, we have 


= [foc SEB eens opr 


where § = 2!-*s—z—e and ¢€ is an arbitrary small positive number and 
S(N) is defined in one of my previous papers.§ 

I shall give elsewhere an application of this theorem to prove that 

G{P(x)} < 17, 
provided that P(x) is a quartic polynomial with coefficient of 2 
positive and that there does not exist an integer g (> 1) such that 
P(x) = P(0) (modq) 

for all x. 

The proof of the theorem depends essentially on the following 
lemma which seems to have some interest in itself. 


+ Proc. London Math. Soc. 41 (1938), 127-42. 
t Ibid. 41 (1938), 143-50. § Ibid. 43 (1937), 161-82. 
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Matn temma. Let P(x) be an integral-valued polynomial of the 


kth degree, and . 
) = >} exp{2mi P(x)a} 
f=} 


1 
Then | | f(o)|- doa = O(pe), 

0 
where {A,u(A)} lies on a polygonal line with vertices (2”,2”—v--e) 
(v = 1.,..., k), and the constants implied by the symbol O depend only 


on the siaihe ients of P(a) and e. 

An improvement of this lemma for the particular case P(x) = x* 
and its application to the additive prime-number theory will appear 
elsewhere later. 

Proof of the main lemma 


1 
Since log( | | f(a)|"d ) 
0 
is a convex function of v,t we only need to prove that 
1 
f(a) |? da = Ofp*—"+*) for v=1, 2...., k. (1) 
0 
Without loss of generality we assume that P(x) is a polynomial 
with integer coefficients. In fact, let g be the least common de- 
nominator of the coefficients of P(x), then, by Hélder’s inequality, 
1 1 a ((p-a)ial . d 
| \f(a)PAda=}|>S >  exp{27i P(qr+a)a}|" da 
0 oe 
1 ((p—a))/¢ 


7 ‘= J z exp[27i{ P(qa~+a)— P(a)} x)!" 


where P(gu-++-a)—P(a) is a polynomial with integer coefficients. 
Then (1) is trivial for vy = 1 and well known for v = 2.8 We are 
going to prove (1) by induction. 
We use the abbreviation 


A Q(x) = 


See, for example, Hardy, Littlewood, and Pélya, Inequalities, § 6.12. 
§ See, for example, Landau, Vorlesungen iiber Zahlentheorie, Bd. 1, Satz 262, 
. There he deals only with the particular case P(x) = a*. For the general 
» see Hua, J. of Chinese Math. Soc. 1 (1936), 23-61, Lemma 11. 
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Then A Q(x) is a polynomial of the (h—1)th degree in x, provided 
that Q(x) is a polynomial of the Ath degree. Let s denote a summa- 
tion with variable x whose number of terms is O(p) 
Consider 
Pp 


Ifa)? = 2 


Z.= 


_ 2 expl2nitP(e,)— P(e,)}a] 
> expl2mi{ Pla. +yi)—P(a)}al 


"i 


Pp 
~ 
= ¥ s exp{2miy, 4 P(aq)a} 


"i ZL 


By Schwarz’s inequality, we have 
Pp p 
|f(«) p> S expl2miy, A Pl, »)at 
Vy x. 


Pp Pp Pp 
<p>>> exp{2miy, ¥2A A P(x 3)a} 

"A V2 Ts 
where A < B means A = O(B). Repeating this process, we obtain 
in general, 

p Pp 4 

f(a) |" << p**-+-1 *.. RS z¢ p(2aiy,...y,A A ...A P(t 41)o} 
Vs Un Ly+4 Yu Vp-1 " 


a ae 3 s * exp{ 2riy; ... 9A -. -AP( Tyusadoy (3 


+p 
"1 Yu Ty+1 


oy ee k—1, where * denotes the condition 


-YyA.. AP(a %y+1) ae 0. 


Vu 


for pu 


We have, therefore, 


1 
f(a) |? da < p?”’-1 J | f(x) |?” * da + 


voy | f(x) |? 


sty i... > > * exp{2miy, .. *Yo-1 yA on! AP 


"” Yy—1 Ly 


0 
By the hypothesis of the induction, the first term on the right of 
(3) 1s O(p”’ ‘1p? *-v+1+e) —_ O(p?’ +€), 
rhe second term on the right of (3) equals 
exp[27i{y, ... Y,— 1, A. -AP(x x,)— 


vee + P(z_>—)}ak] da = pp 


Dp 


pi... 
2' »R, 
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where R denotes the number of solutions of 
Y,---Y,-1 A A... P(z,) = P(z,)—P(z-)+...—P (2), 


Yui Via 
Y *Yv-1, A A P(x y) # 0, eu > Uw * ey 41 <P. (4) 
For given 2,,...,2»-1, the number of solutions of (4) is 
Ofd’-"[ P(z,)— P(z.) +... — P(z_) ]}.T 


Since d(n) O(n£), we have 
R.: > cee > t@&-{P(z,) —_ P(2,)+... — P(z,» 1)} 


p ad 
where { denotes the condition P(z,)— P(z,)+...—P(z»—) # 0. The 
lemma is therefore proved. 


Proof of the theorem 
Let p = N* and 
Pr 
Sia) = > exp{27i P(x)a}, 
zr=1 
where p, is the greatest root of P(x) = N. It is evident that when 
N is sufficiently large, p, — exists and p< p,<p. Then 


‘| IIs xjexp(—27ialN) da. 


It is sufficient to waieaan the part W of the integral corresponding 
to the minor arcs, since the remaining part can be treated, without 
any difficulty, by the method used in my previous paper.§ 

By Weyl’s theorem!| and the main lemma we have 


W <pi-2™ reXs—2" {ts (x) )| dex 


ak 1 ° 
g1-*, eyg_ga*) i ok 
pa-2*+ens-29( | [ [ |S(c)|2 dx) 
r . 
0 


pk, 


The theorem is proved. 
In closing, I should like to express my warmest thanks to the 
referee for his valuable advice. 
+ d(n) denotes number of divisors of n. 


§ Proc. London Math. Soc. 43 (1937), 161-82. 
Landau, J adiiaden liber Zahlentheorie, Satz 267. 











A PAIR OF FUNCTIONS WHICH ARE FOURIER 
SINE-TRANSFORMS OF EACH OTHER 


By R. 8. VARMA (Cawnpore) 
[Received 22 February 1938] 
Tue object of this paper is to establish the following 


THEOREM. If 
f(x) = KD® 6 (02) —D® n.4y(—12)} 


9 
- Wa) = (—" Fe L, (a), 
nN: 


then (x) and x(x) are Fourier sine-transforms of each other provided 
that n is a positive integer. 

The following two lemmas are required: 

Lemma 1. If n is a positive integer and m > —}, then 


= | 5; = ea 
n u l Tm+Uq) 4 1.3.(n+1)(n-+-2) [+2 n)— 


m I 9 4 m 


Tn (2) =, 


m 
In a recent paper* I have shown that 
re n 
e brgmJ™ (2) : p(s Pp) 
’ m\" n! (4+ pment 


(1 —p)" 


This gives a" T%, (x) = 5 
n y™ +n 


m 


by the help of the operational formula, 


e-*f(z) = P d( p+), 


pra 
where f(x) = o(p), 
Now consider the series 


(1—p)"{,__ 41-9 (=). \_ 1 (1—p)" 


! m+n | ae +H ‘| aot eee 
n! p i! p 2! \ p n! p 


If we interpret it term by term, we at once obtain the above lemma. 


* R. S. Varma, ‘Some functions which are self-reciprocal in the Hankel 
transform’: Proc. London Math. Soc. (2), 42 (1937), 9-17. 
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9n+h 
Since Th ) oe (2n-+-1)! yet te! Dean saly (2)], 
this gives, for m = }, 


Wm) ayia 
(2n+1)! Pon (anv 2)—5 


1.3 
+ ~ an+5(UV 2 
2.4.(2n+3)(2n+5) 7"* 


Lemma 2. If n is a positive integer, 


(—)™i{D? (n+l ix|(22a)]—D? (, epL—ty (2ar)]} 
gn+l \/(Qar ){ 1 
= P. Vx D,,, +9(2Vx 
(2n | 1)! | 2n +3 (2 )+ | 2(2n +3) 2n 3( )+ 
‘ de Don 4s(2V2)+.. 
" 2.4.(2n+-3)(2n+5) *”*” 
rhis has been proved by 8. C. Mitra.* 
lo establish the above theorem 
even integral values of n. 


7 


, let us first restrict ourselves to 
We then have, using (2) 


8 
a 


9) 
[(-) | 2re-V T"(y?) sinay dy 
7 
7 0 


+1 /(27r) 
—¥ )! '{\p 2n+ ove ee Dz, .g(y¥V2) + 


Don 45(YV2)— tein xy dy. 
) J 
Integrating term by tern 


process obviously justifiable—and 
making use of the formulat 


nN Pe 
(") | Domsa(y¥2)sin xy dy = (—)"Day(02), 
N \7/ . 
0 


* §. C. Mitra, ‘The operational representation of D,(x) and 
{D* (n41)(tx) —D* (ny)(—ix)}, 
Proc. Edinburgh Math. Soc. 4 (1934), 33-5. 
+ A. Milne, ‘On the equation of the parabolic cylinder function’: Proc. 
Edinburgh Math. Soc. (1), 32 (1914), 2-14. See also E. C. Titchmarsh, Intro- 
duction to the theory of Fourier Integrals (Oxford, 1937), 261 
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we obtain 


9) > 
J) | 2re-*¥* 7" (y?)sin zy dy 
7 
0 


2n+1, (Qa) | l 
an tol ) 
Gert | ansil®82) +9 (on 13) 
1.3 
+ 2.4.(2n+-3)(2n-++5) 
= UD? (nn (t2)—D? (yn 4y(—t2)} 


= $(2), (3) 


where we have used Lemma 2. 


Dy», g(X v2) + 


Dayss(t2)+--| 


Proceeding in the above manner, it can be shown that 
2\ - . 
7 | d(y)sin vy dy = 2ne*” T"(x?). (4) 
NWT 
0 


It is also easy to see that the results (3) and (4) remain true when 
nm is an odd integer. Moreover, Sonine’s polynomial is related to 
ordinary Laguerre polynomial by means of the relation 

ies (— l * 

T 3 (x) — L,, (x). 


n! 


This proves the theorem stated above. 











SOME THEOREMS ON FUNCTIONS REGULAR 
IN AN ANGLE 
By V. GANAPATHY IYER (Madras) 
[Received 27 April 1938] 

1. I HAVE shown* recently that there exist sequences of points which 
resemble closely the lattice-points 

z=m+in (m,n = 0,+1, +2....), (1) 
so far, at any rate, as the behaviour of integral functions at these 
points are concerned. Many writers} have considered the properties 
of functions regular in an angle in relation to the lattice-points in 
that angle. The main object of this paper is to examine whether the 
allied sequences mentioned above possess similar properties in rela- 
tion to functions regular in an angle. 

1.1. Let z = re’ and let A(a) denote the region 
Ol\<a (0<a<n), (2) 
where « is fixed in any particular discussion. Let f(z) be a function 
regular in A(«) except perhaps at z = 0, where it is supposed to be 
continuous as z +> 0 from within A(a«). We write 
M{r,f, A(x)} = max|f(z) 

as z varies in the sector common to A(a) and the circle |z| < r. The 
order p and the type «{f,A(«)} of f(z) in A(a«) are defined by the 


relations 


loglog M{r, f, A()} 
ee log r ; ra rP 

When f(z) is an integral function we write M(r,f) = M{r,f, A()} and 

«(f) = «{f,A(z)}. We shall denote by C(p,d) the class of all integral 

functions of order p and type less than d, that is, all those functions 


: log Mfr, f, A(«)} 


p = lin K{f, A(x)} = lin 


f(z) satisfying the condition 
—log M(r,f 
lim 2 (7,9) < d, 
r—> oo rP 

where p and d are two finite positive numbers. 

* Trans. American Math. Soc. 42 (1937), 358-65. This will be referred to as 
E in the sequel. 

+ M. L. Cartwright, Proc. London Math. Soc. (2), 43 (1937), 26-32. A. 
Pfluger, ibid. (2), 42 (1937), 305-15. J. M. Whittaker, Proc. Edinburgh Math. 
Soc. 2 (1930), 111-28; Proc. London Math. Soc. (2), 37 (1934), 383-401 ; Inter- 
polatory Function Theory (Cambridge Math. Tracts, 33, 1935). 
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1.2. We now proceed to specify the kind of sequences to be con- 
sidered in this paper. Let [z,] (z, = 7,,e") be a sequence of distinct 


complex numbers such that 
0<', <%=. =, 70 


as n->0o. For shortness we shall refer to the exponent of conver- 
gence of [z,] as its order. Let B be a positive constant and let 
A,,(h) denote the circle with centre z, and radius Bj|z, |-”. We denote 
by A(h) the system of circles {A,(h)}. Let p and d be two finite 
positive numbers. We shall say that the sequence [z,| belongs to 
the class L(p,d) when there exists an integral function g(z) of order 
p with simple zeros at the points of [z, | (and having no other zeros) 


satisfying the two following conditions: 


(i) lim logig’@n)| = d, 
nx Zn |? 
... log|g(z)| (3) 
(ii) = >d 
as |z| > 00 outside the circles A(h) for some positive B and some h 
greater than p. 

1.3. The lattice-points (1) are the simplest known set of the type 
discussed in § 1.2. It belongs to the class L(2, $7) and (ii) of (3) holds 
for any positive B and any h greater than 2. It is easily shown by 
using classical theorems* on integral functions that any sequence 
[z,,| for which (ii) of (3) holds is of order p so that the sum of the 
diameters of the circles of A(h) is finite. I have proved (E, Theorem 1) 
that, if [z,,] belongs to L(p,d), the set [8,,] of the amplitudes of the 
points of [z, | is everywhere dense in 0 < @ < 27, so that every angle 
however small contains an infinity of points of [z,]. The sequence 


28.0 (4) 


belongs to L(p,d). The method used to prove that (4) belongs to 
L(p,d) can be easily modified to prove the following lemma which 
gives a general class of sequences belonging to L/(p, d). 


(E, Lemma 4) 


| (2pd)- Upp 2ipg2vmi)in} 


* Cf. my papers in J. Indian Math. Soc. (New Series), 2 (1936), 4-6 and 
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Lemma 1. Let I(x) and m(x) be two differentiable functions of x in 
0 <2 < ow tending to zeroas x > 0. As x >, let 
al'(x)>0 and xlogxm'(x)>0. 
Let «, n, A, w, L, and M be positive constants. Let 
[A,] (O<A, <A... <A, > 0) 
be a sequence defined by 
A, = ImM1+l(n)}. (5) 
Let [p,| (0 < py < po < pg.» <p, > 0) be a sequence of integers 
defined by By, = Mne{1+m(n)}. 
Then the sequence 


| nas e2v7D/pn| ( 


belongs to L(p,d) where 


n 


V 


_#t1 

Oe ee mP(e + 1) Dr 

1.4. Given p, d, L, M, A, », the numbers « and 7 can be so 

determined from the relations (8) that (7) belongs to L(p,d). The 

main properties of integral functions in relation to a sequence of 
L(p,d) are contained in the following theorem. 

THEOREM 1.* Let f(z) be an integral function of the class C(p, d) 
[$1.1] and let [z,,| be a sequence of L(p,d). Then 


lim 


n—- DB 


«(f) = lim 


TD 


log M(r,f) 
rP 


log |f(2,,) 
~ p F 
! (7 


‘ 


Also, a function of C(p,d) bounded at the points of a sequence of L(p, d) 
must be a constant. 

1.5. In this paper I propose to prove results similar to those of 
Theorem 1 for functions regular in an angle. The method used 
closely resembles the one I have adopted elsewheref to derive the 
asymptotic behaviour of a function along a straight line when its 
behaviour at a sequence of points on the line is known. The special 
methods that have been used in the case of the lattice-points by the 
previous writers cannot be applied here, since these make use of the 
peculiar properties of lattice-points considered as complex integers. 


* KE, Theorem 2 and footnote on p. 365. 
+ Quart. J. of Math. (Oxford), 8 (1937), 131-41, Theorem 1. 
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In Theorems 2 and 3 below, the method used here has necessitated 
the assumption that the sides of the angle A(a) have no point in 
common with the circles of A(h) except perhaps with a finite number 
of them. It is very likely that this assumption is quite unnecessary, 
since it does not appear to have any bearing on the remaining hypo- 
theses and conclusions of these theorems. I show later (Theorem 4) 
that in the case of the lattice-points (1) and the sequences (4) and 
(7) this hypothesis can be dispensed with; in fact, I show that this 
is so in the case of a definite sub-class of L(p,d). But I have not been 
able to get rid of this assumption in the general case of a sequence 
of L(p, d). 

2. The two following theorems constitute the main results of this 
paper. 

THEOREM 2. Let f(z) be regular in A(x) and let [z,,| be a sequence of 
I(p,d). Let A(h) be the system of circles in (ii) of (3) and let the following 
conditions hold: 

(i) the sides 0 = +« of A(a) have no point in common with A(h) 
except perhaps with a finite number of these circles; 


log M{r, f, A(x)} 
rP 


(ii) lim <8 


r—>o 


(iii) lim 
T—> 00 


log |f(re***)| _ B <d; 
rP — > 


x | | 
(iv) lim log fen) <B<d {z,cA(a)}. 
cilities =n | 


Then {fb (a)} = fim 8M ey A} < g <a, (9) 
THEOREM 3. Let f(z) and [z,] be as in Theorem 2. Let the conditions 
(i) and (ii) of Theorem 2 hold. Let 
(iii)’ |f(re***)| < B; 
(iv)’ |f@,)| < By, for 2, cA(q). 
Then, for all z in A(x), If(z)| < B. (10) 
2.1. We shall first establish a few lemmas. Lemma 2 has been 
proved in E, Lemma 2. 
Lemma 2. The function g(z) in the definition of the class L(p,d) 
[§ 1.2 (3)] as of order p and type d. 


3695.9 Pp 
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Lemma 3. Let A > 0. Then any ring R—A < |z| < R contains a 
circle not cutting any circle of the system A(h) provided that h > p and 
R>R, = R,(). 


Proof. This follows from the fact that > |z,|-* converges when 
h > p. 


Lemma 4. Let the sides of A(x) have no point in common with A(h). 
Let the remaining conditions of Theorem 2 hold. Let B < y << dand 


f(z) f@n) 1 
H z) = zm -_ gm : ll 
on ) @ g(2) n g' (Zp) z—2 ( ) 


Zn A(x) - 


where m > 0 is an integer and g(z) is the function in the relation (3). 
Then, for all z in A(«), 


|Hn(2)| 


< AU. AO. 


mc A(x) 


ml|p 
| |jm+hexpl —(d— le |p rs 
12m | exp[ (d y)| n| +s | |, 


(12) 
where A(f,y) is a constant* depending on the function f(z) and the 
number y. 


Proof. Since 6 = +« have no point in common with A(h), we have, 
for z on 0 = +a, 
tnd(Z,) 1 


te poe | SAGYAGY) DY lenl™**exp[—(—y)l2nl?], 
g (2) ov 2nc A(a) (13) 


in virtue of (i) of (3) and the condition (iv) of Theorem 2. By (ii) of 
(3) and the condition (iii) of Theorem 2 we get, for z on 6 = +a, 


amd (2) 
g(z) 


| Sync A(x) 


| <A(f, y)A(g, y){ max r™ exp[—(d— yr?) 


< Alf, yAG,y)|—™ ¥ 
Un »| a 


* Throughout the rest of this paper an expression of the form 


A(f,g a, 
will denote a finite positive constant depending on the functions f, g,... and the 
numbers y, A,... which need not have the same value at the different places 
where it occurs. It is supposed that the number B occurring in the radius of 
the circles of A(h) is unity ; otherwise B~' will occur as factor in (12) and this 
does not affect the arguments in § 2.2. 
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By (ii) of Theorem 2 there is a 8 (0 < 8 < d) and a sequence 
[R,] (0< R, < R,... < R, >), 
such that M{R,,, f, A(«)} < exp[(d—8) Re]. (15) 
By Lemma 2 there is a circle in the ring R,—1 < |z| < R, which 
does not cut any circle of A(h). If Tj, denotes ay are of this circle 
— \f(@)| < expl(d—3) Rf] (16) 


for z on [, in virtue of (15). Since [, is outside the circles A(h), we 


can use (ii) of (3) and in conjunction with (16) conclude that 

gm g(2) 

f(2) 

The relations (13), (14), (16’) give (12) in conjunction with the 
maximum-modulus principle. 


+0 as |z|—-oo along the ares I. (16’) 


2.2. Proof of Theorem 2. If we omit a finite number of points from 
a sequence of L(p, d), the rest still belongs to L(p,d). Therefore we can 
suppose without loss of generality in the condition (i) of Theorem 2 
that no circle of A(h) cuts @ = +a. We can now apply Lemma 4, 
and so the relation (12) is valid for any y such that B < y<d. 
Let 0 < 8 < d—f, and choose A so that dit = d—B—8 (if B = —oo 
we can choose A to be any large positive number). Let 


x(z) = Co +e, 2+-€,27-+-... 
be any integral function of order p and type not exceeding d)/. 


Then, by a classical theorem on integral functions, 


[Nm < (dArep)*?. (17) 


lim m¥P|c 


Take y = B+45 = d—ddP—}6 in (12). Then we get in virtue of (i) 
of (3) and (17) 
ee flen)x(@) 1 
g' (Zn) Z—2y 


m 


onc A(x) 


<A VAG YE lenl*expl—(P+18) 2 °1( ¥ ml én) 


se ; esl 
+ > eden 
y)A(x, p, d, A, 5, A) 


? iii. Pog ° (18) 
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By Lemma 2 we can take x(z) = g(Az) in (18). Then we get for z in 
A(x) but outside the circles A(h), 
| (2) | gop 
z)|< A(f,g,p,d, B, 5, h), 19 
f(2)| g(Xz) (f,9,p,4,B, 5, h) (19) 
in virtue of (i) of (3) and the choice of A. Let A)(h) denote the system 
of circles round the zeros of g(Az) similar to A(h). Then by (ii) of (3) 
and (19) we get that 
—log|f(z)| 
lim es < d—dn = B+45, (20) 
as |z| > 00 outside the circles of A(h) and A)(h). Let « > 0. We get 
from (20) that for any z outside these circles 
f(2)| < exp[(B+-46+)|z|?] (21) 
provided that |z| > 7) = ro(e). Now the system of circles A(h) and 
A,(h) together determine a sequence of non-overlapping domains 
D,, Dg,..., Ds,..... The diameter of any D, does not exceed the finite 
number P=2 > Biz,|-(1+a~). 
nc A(x) 
Since (21) holds on the boundary of D,, for all nm > no, and since the 
diameter of D, does not exceed P, it follows from the maximum- 
modulus principle that, for z inside any D.,, 
If(2)| < exp[(B-+ 38426) |2/P] 
provided n > n, = m,(e). From (21) and (22) we conclude that 
log | f(z)| 
olf)! < p+ 4s, (23) 


lim 


as |z| > 00 in A(a). Since 6 is as small as we please in (23), we get the 
relation (9). This proves Theorem 2. 

2.3. Proof of Theorem 3. As in the proof of Theorem 2 we suppose 
that 6 = +a have no point in common with A(h). Theorem 2 in 


conjunction with (iii)’ and (iv)’ of Theorem 3 gives 


y Se fF / 
, log Mf) SMe} < 0. (24) 
r 


«ff, A(o)} = lin 


Tro 
Hence (24) holds when f(z) is replaced by | f(z)|”, where p > 0 is any 
integer. Proceeding as in the proof of Lemma 4 and using (iii)’ and 
(iv)’ of the present theorem we get, for all z in A(a), 


Lie? _ > Lied 1 | < 4G,a)Be+AG,d, WBE. (25) 
g(z) J (2) Z—Zy 


Snc A(x) 
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From (25) we get, for z in A(a«) but outside A(h), 
p 
Ue < A(g,d)B”+2A(g,d,h) Bp, . 
from which we get, for any fixed z in A(«) but outside A(h), 
\f(2)| < lim |g(z)|"”[A(g, d)B” + 2A (g, d, h) Bp}? 
pa 
< max(B, B,). (26) 
Using the maximum-modulus principle as in the last stages of the 
proof of Theorem 2, we find that (26) holds throughout A(«). Hence 
f(z) is bounded in A(a) while |f(re*'*)| < B, so that by applying a 
well-known theorem of Phragmén and Lindeléf* we find that (10) 
holds for all z in A(q). 

3. We shall now show that in certain cases the condition (i) of 
Theorems 2 and 3 can be dispensed with. For this purpose we 
require a preliminary result. The circles A(h) in (ii) of (3) determine 
a sequence of non-overlapping domains 

[ D,, D,,..., D,...] = [D,] = Dia). 
Let D,(h) be a sub-sequence of domains selected from those of D(h). 
We shall prove the following 

Lemma 5. Let [z,] be a sequence of L(p,d). Let those points of 
[z,, | Lying in D,(h) be a sequence of order less than p. Then the sequence 
of points of {z,,| not lying in D,(h) belongs to L(p, d). 

Proof. Let the points of [z,] lying in D,(h) be of order py < p and 
let o(z) be the canonical product with simple zeros at these points. 
Then, if py < po < p, we havet 

log|o(z)| = O(|2\P*), (27) 
as |z|—>0o outside D,(h). Let D,(h) denote those domains of D(h) 
not belonging to D,(h) and let 


G(z) = 9) 


oi o(z) . 


Then G(z) is an integral function having simple zeros at those 


* G. Valiron, Lectures on Integral Functions, p. 125. 
+ Here we use the fact that, if g(z) be a function of order p and A(h) denote 
the circles round the zeros of g(z) as in § 1.2, we have 
log|g(z)| = O(|z|") (p’ > p) 
as |z|—> 00 outside A(h) for any finite h. This result can be obtained by argu- 
ments similar to those used in the usual proof of Hadamard’s theorem on 
the minimum modulus of integral functions. 
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points of [z,,] lying in D,(h) and having no other zeros. In virtue of 
(27) and the fact that the domains of D(h) are non-overlapping we 
get, for z, c D,(h), 
log |G" (z,)| = log|g’(z,)|+-O(|2n|"); (2 

so that (i) of (3) holds for the points of [z,] in D,(h). Also, by (ii) of 
(3) and the fact that D,(h) c D(h), we get 

log |@(z)| = log|g(z)|+-O(|z |?) (29) 
for z outside A(h) = D(h). Hence 

|G(z 
lim 2814 a d, (30) 

as |z| oo outside A(h) = D(h). Since G(z) does not vanish in 
D,(h), log|G(z)| is harmonic therein. Using the fact that a function 
harmonic in a domain D cannot attain its maximum or minimum in 
the interior of D and arguing as in the last stages of the proof of 
Theorem 2 we can conclude that (30) holds as |z| > 00 in D,(h) as well. 
Hence (30) holds as |z| > 00 outside D,(h). This implies that (ii) of (3) 
holds for points of [z, | lying in D,(h). This proves the lemma. 


3.1. Now let [z,| be a sequence of L(p,d). Let A(h), D(h) be as 
in §3. Let [z,,4, D(h)| denote those points of [z,,] that lie in such 


domains of D(h) as have a point in common with the line @ = 4). 
If for every 0 in 0 < 6 < 2m the sequence |z,,0, D(h)| is of order less 
than p, we shall say that |z,| belongs to the class L,(p,d). It is easy to 
see that the lattice-points (1) and the sequences (4) and (7) belong to 
L,(p,d). Using Lemma 5 we see that a sequence [z,,] of L,(p,d) still 
belongs to the class when the points of [z,,, «, D(h)| and [z,,, —a, D(h)] 
are omitted from it. Hence in the case of such sequences the condi- 
tions (i) of Theorems 2 and 3 can be omitted and the theorems still 
remain true. So we can state 

THEOREM 4. Let [z,] be a sequence of L,(p,d). Then Theorems 2 
and 3 hold even when conditions (i) of these theorems are dispensed with. 
In particular this is true in the case of the sequences (1), (4), (7). 

3.2. It is a significant feature of the results of this paper that the 
magnitude of the angle a of A(«) does not affect the truth of the 
theorems proved here. If a < 7/2, the conditions (ii) and (iii) of 
Theorem 2 and (ii) and (iii)’ of Theorem 3 are sufficient to ensure 
that (9) and (10) hold in the respective cases. But, if «a > 7/2p, this 
is no longer true and in this case the conditions (iv) and (iv)’ come 
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_ into play. In this connexion it is interesting to compare the follow- 
ing result due to Miss Cartwright.* 


THEOREM 5. Let f(z) be regular in A(«). Let the following conditions 
hold: 
id 
(i) lim —_ U <A< }r for a—8 < |6| < a for some positive 5 


To 
and X < 4rsin®a tf « < 4; 
(ii) lim log Mir. f, A(a)} <hr; 
ra r 
(iii) f(m+-in) = O(1) for all those lattice-points lying in A(a). Let 
B < « be such that sin®(a—B) > 2A/7. Then f(z) = O(1) for z in A(B). 
3.3. The above theorem states that under the conditions (i), (ii), 
and (iii) of the theorem there are angles inside A(a) in which f(z) is 
bounded. Since the hypotheses are with reference to A(«) and the 
result proved is for an angle A(8) (8 < «), a condition similar to (i) 
of Theorem 3 is fulfilled effectively. In Theorem 3 condition (iii)’ 
replaces condition (i) of Theorem 5, which is less stringent. But, as 
shown in Theorem 4, the condition (i) of Theorem 3 on the boundary 
of A(«) is unnecessary in the case of lattice-points. It is very prob- 
able that a result of the type in Theorem 5 holds for any sequence 
of L(p,d), though the methods of this paper seem insufficient to 
prove it. But Theorem 2 of this paper enables us to replace (i) of 
Theorem 5 by the less restrictive condition 


| +ia)| 
lim logifire Pei< te (31) 
ro ad 
since (31) in conjunction with (ii) and (iii) of Theorem 5 implies that 
fe £ Al, 
— log Mir. J, A(a)} <i<ch 
ro r 


and the last relation contains, a fortiori, the condition (i) of Theorem 5. 


[Added 21 July 1938] 

A step has been omitted in the proof of Lemma 4. In addition to 
(16’), it is necessary to show that the expression on the left of (13) 
tends to zero as |z| > 00 along the ares [},. But it is easy to see that 
this expression tends to zero as |z| > 00 outside A(h) by considering 
the two sums corresponding to |z,,| < (1—A)|z| and |z,| > (1—A)|2| 
where 0 <A <1. A fortiori the same is true as |z| > 0 along [,. 


* See M. L. Cartwright, loc. cit. 








ON DIVISOR PROBLEMS 
By E. C. TITCHMARSH (Ozford) 
[Received 26 January 1938] 
1. Let d,(n) denote the number of expressions of n as a product of 


k factors, and let D,(a) = ¥ d,(n). 


n 


It is known that 


D(x) = (Qo 4+41% log x+-...+-a,,_, ;, log*-tx)a+-A,(x) 

where A,(2) = o(a). We may define the order a; of A;,(x) as the 

lower bound of numbers A such that A,(#) = O(x*). The problem of 

determining a, is notoriously difficult.* It was proved by Hardyt+ 

that a, > 3(k—1)/k, but the exact value of a; is not known for any k. 
An easier problem is that of determining f;,, the average order of 

A,.(«), defined as the lower bound of numbers A such that 


: | Az(y) dy = O(x?). 
x 
1 


Clearly 8), < a. It was proved by Hardy? that 8, < 4; by Cramert 
that 8, < 1—2/k for k > 3; and by Hardy and Littlewood§ that, 
on the unproved Lindeléf hypothesis, 


By = H(k—1)/k 
for every k. 
All these results depend more or less on the Riemann zeta-function. 
The connexion is given by the formula 
{C(s)}* = M3 =e io Zp 
n=1 
Apparently no one has pointed out explicitly the exact relation 
between the f,, and the theory of the zeta-function. It is, however, 
quite simple. Let o; be the lower bound of values of o for which 
r 
| |S(o+-it)|%* dt = O(7). 
1 
Then, for each integer k, a necessary and sufficient condition that 
By = 3(k—1)/k ts that o, < ¥(k+1)/k. 
* See Hardy and Littlewood (5). t+ Hardy (3), (4). 
{ Cramer (1). § Hardy and Littlewood (6). 
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I also prove that Bf, = }, 8, =}, that ?< 6, <3, and that 
B;, > 4(k—1)/k for every k. 


2. Perron’s formula for the sum of the coefficients in a Dirichlet 


series gives 
1 e+iT m 
D,{x) = —. lim | ¢ (8) 16 ds (c>1). (2.1) 
271 To 8 
e—tT 
Applying Cauchy’s theorem to the rectangle y—iT7’, c—iT, c+iT, 
y+iT’, where y is less than but sufficiently near to 1, and allowing 
for the residue at s = 1, we obtain 


(2.2) 


Let y;, be the lower bound of positive numbers o for which 


\f(o+-it) |?* 
io tit? -dt <<. .3) 


. 


for o > 0, it follows from a general theorem on mean-values of 
analytic functions* that the left-hand side of (2.3) is a bounded 
function of o for y,+e <o<1—e (e > 0). Hencet C*(s)/s +0 
uniformly as ¢- -oo inside any strip y;,+« < o< 1—e, where 
« > 0. Hence, if we integrate the integrand of (2.2) round the 
rectangle y’—iT’, y—iT, y+-1T, y’'+iT, where y, < y’ < y <1, and 
make 7' > 00, we obtain the same result with y’ instead of y. 

If we replace x by 1/x, (2.2) expresses the relation between the 


Then (2.2) holds for y, << y <1. Since {(s) = O(|t|4) uniformly 


Mellin transforms 
_ has Xie) — Lk(g 
fix) =A,(/z), —&(s) = L*(s)/s, 
the relevant integrals holding also in the mean-square sense. Hence 
Parseval’s formula for Mellin transformst gives 
ow 


a “4) | 2k 
Lf er t= [ ae(;}e dec = | Ap(x)u-2°-1 dar, (2.4) 
2a jo-+-at|* J xv 

0 


. 
— 0 


. 


provided that y, <0 < 1. 


* Hardy, Ingham, and Pélya (7), Theorem 7. 
+ See, e.g., Titchmarsh (9), § 5.4, lemma. 
t See Titchmarsh (9), Theorem 71. 
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It follows that, if y, << o < 1, 

x 

| Az(x)a-2e-1 dx < C = C(k,o), 

4X 
x 
AR(x) dx < CX2e+1, 
x 

and, replacing X by 3X, 1X,... and adding, 

x 
| Az(x) da < CXe+1, 
i 


Hence f, < o for all such o, i.e. By, < yx. 
The inverse Mellin formula is 
a 


a(t) tdz = [ A,.(x)a-$-1 dx. (2.5) 


« 


0 


The right-hand side exists primarily in the mean-square sense, for 
VY, <o<1. But actually the right-hand side is uniformly con- 
vergent in any region interior to the strip 8, << o <1. For 
x 


x x 3 
20-2 dy} * 
J 


| 
tx xX pX 


— {O(X *Br+1+€)O(X-2e-1)\3 res O( XBi-o+t€) 


2 x J 
( |A,.(a) |a-°-4 da < [ A}(x) dx | 


and on putting X = 2, 4, 8,... and adding, we obtain 
| |A,(x)|a-°-2da<C (B.<oa< 1). 
i 
It follows that the right-hand side of (2.5) represents an analytic 
function, regular for 8B, << o«< 1. The formula therefore holds by 
analytic continuation throughout this strip. Also, by the argument 
just given, the right-hand side of (2.4) is finite for B, <o< 1. 
Hence, by the Plancherel theory of Mellin transforms, so is the left- 
hand side, and the formula holds. Hence y, < f;, and so in fact 


Yr = By 
3. We can now prove that B, > 4(k—1)/k. If } < o < 1 we have 


P T u(2k)/ FT \1-1/(2k) 
CT < | \f(o+it) |? dt < ( \f(a-+-it) |* i) | a F 
iT iT 


iT 
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T 
whence [ \f(o-+it)|** dt > C’T. 
1T 
Hence, if0<o0<4,7>1, 


T 
r If(o-+it)| . S(o+i|* 4 C a 
jo--it|? u> [i lob it |? t> a | [S(o-+at) | dt 


iT 
> orpea.sn-s f |f(1—o—it) |** dt > CT*1-20)-1, 
ir 
This can be made as large as we —_ by choice of T if o < }(k—1)/k. 
Hence y; > $(k—1)/k, ie. B, > 4(k—1)/k. 
To prove the theorem stated in § 1, suppose first that o, < $(k+-1)/k. 
Then, by the functional equation, 
f | j 
| |(o-+-it) |* dt = ss aan |(1—o—it) |* «) == Q(FU-sn) 
i i 
for o < 4(k—1)/k. It follows from the convexity of mean values of 
analytic functions* that for 4(k—1—e)/k < o < }(k+1+6)/k 
yr 
[ g(o-+in) 2 at = O(7 
i 
This is of the form O(T?-*) (8 > 0), if o > 3(k—1+e)/k. Then 


1+ (3+ T 9) 


- 
| \f(a-+-it) |? dt = O(T-), 
J |jo+-at|? 


4T 
and it follows that the left-hand side of (2.4) is finite. Hence 


ve < Hk—Ifk. 


Hence y, = B, = 3(k—1)/k. 
On the other hand, if 8, = $(k—1)/k, it follows from (2.4) that 
fi £(o+it)|** dt = O(T?) 
for o > $(k—1)/k. ecg! by the functional equation, 
fi IC(o-bit) | dt = O( T1-20)-2) 


1 


* Hardy, Ingham, and Pélya (7), Theorem 6. 
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for o < }(k+1)/k. Hence, by the convexity theorem, the left-hand 


side is O(7!+*) for o = $(k+1)/k, and so, by an argument due to 
Carlson,* is O(7') for o > $(k+1)/k. Hence o, < $(k+1)/k. 

4. It is known} that o, < 1—1/k. Since 1—1/k < }(k+1)/k for 
k < 3, it follows that B, = }, B, = }. 

The available material is not quite sufficient to determine £,. We 
have 8, > 3. To obtain an upper bound for it, we observe that, since 
C(3+it) = O(t**), 

. ~ 
[ \€(4+- at) |§ dt = o(ri+« [ icin ita = O(T#**), 
1 1 


Also, since o, < %.t 
1 


r 
\C(3,+7t)|§ dt = (1% 


1 i 


C(5—7)|§ ay — O( T3/5+€), 


Hence by the convexity theorem 
r 
| [S(o+-it)|§ dt = O(T*-¥e+*) (<a < }). 
1 
This is of the form O(T?-*) if ¢ > 3, and it follows as before that 
B, <3. A slightly better result could be obtained by using more 


refined results about ¢($+-7?). 


* See Titchmarsh (8), § 2.42. + Titchmarsh (8), § 2.51. 
t Davenport (2). 
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ELECTRICAL NOTES 


By F. B. PIDDUCK (Ozford) 
XI. MAGNETRON OSCILLATIONS* 
[Received 11 March 1938] 


12. Calculation of current for the first orbit 
The convection current due to a single charge e is de,/edt, from 
(14), where (apart from a constant) 


7e,/e = imaginary part of [—log(i—z/b)+log(—i—z/b)] (36) 
and z/b = em}(wt—at,)exp i[ 6, + }(wt—wt,)]. (37) 
The charge grazing the anode, in the steady motion, between polar 
angles 0, and 0,+d6, and between times ¢, and ¢,+dt, is 


—clI d6, dt,/27; 
and, as in Note VIII,+ the current due to all the charges emitted so 
as to graze the anode between 6, and 6,+d@, is 
I d0,(f.—f,t+Sa—fet ---)/27, (38) 
where f stands for the right-hand side of (36) and the electrons 
between the electrodes, beginning with those which left first, extend 
at time ¢t from 2, to z,, from z, to z,, and so on. 

Let the period of the oscillations be the same as the time of one 
orbit, let b/a = 500 or v/w = 0-903 as before, and let the range of 
vt, for anode capture be from a to 7+, where a is known for every 
value of 6,. Remembering that an electron caught by the anode is 
present in the field when —z < vt—vt, < 0 and that one not so 
caught is present when —z < vt—vt, < 7, it is easy to see that 
electrons are present when vt < vt, < 27+-a if a < vt << 7+a, and 
when 7+a < vt, << 7+ut if m+a< vt < 27+a. That is, there are 
only two terms in equation (37), such that 


Re vt— oa P vt — 
2 = om exp él 6,+ 


"1 _ exni,. a 
"iis ili b 1-806 1-806 


when a < vt < a+a, and 


> a ek 
1-806 b 


b 1-806 


* This note is a continuation of Notes IX and X, Quart. J. of Math. 
(Oxford), 7 (1936), 241, and 8 (1937), 75. 
+ Quart. J. of Math. (Oxford), 7 (1936), 211. 


“— exp i + 
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when 7+a < vt < 27+a. Let 0, lie at first between 0 and 7. Then 
when a < vt < 2+a, f, is $ or 3 according as 6, < or > 41. If 0, is 
increased by z and vt is unaltered, f,—/, is replaced by f,—1 when 
x<vt<a+a and by 4$—f, or 3—f, when a7+a < vt < 27+, 
according as 6, < or > 4m. The calculations, therefore, need only 
be carried out in full for 0 < 6, < 7; and as 0, = 47 has to serve for 
a range of values round $7, the fractions $ and 3 are then replaced 
by 1. 

The table on pp. 223-4 gives the values of 1000(f,—/,) for 6, = 0° 
to 350° and vt = 0° to 170° at intervals of 10°, the values repeating 
with reversed sign when 6, and vt are increased by 180°. This table, 
therefore, gives the ratio of the convection current to the emission 
current (multiplied by 1000) for a fine stream emitted from the fila- 
ment at any assigned angle. The present method is more general 
than that of § 5, and can be applied to undisturbed orbits of any 
form when electrons are emitted at a constant rate from all the sur- 


face elements of the filament, provided that they are all caught when 
they return. Treating the columns of the table by Simpson’s rule, 
we have the corresponding ratio for a filament emitting uniformly; 


numbers which are given in the last row on p. 224, below the main 
entries. These currents differ considerably from those of § 5. The 
difference is due to three main causes. Firstly, there is the approxi- 
mation of § 4, where the range of 6, for capture by the anode was 
taken to be 180° for all values of wt,. To see what error it introduces, 
I calculated the entries under wt = 0° and 90° in the table of § 5 by 
the new method and found —0-209 and 0-103 instead of —0-222 and 
0-096. The error is thus of the order of 7 per cent. The remaining 
difference is due partly to precession and partly to the strong field 
near the edges of the sectors, which alters the limits of vt, for capture 
by the anode, especially near 0, = 90°. A magnetron with a narrow 
gap will differ appreciably from one with a wide gap. Since the gap 
is often quite wide, it may actually be nearer to the truth to suppose 
that the disturbing field is uniform, though no accurate theory is 
possible for a wide gap.- The fundamental term in the convection 
current taken from the last table is the real part of 


(—0-204—0-1917)/ exp wt, (39) 


instead of (—0-233—0-0987)/ exp iwt in equation (15), The coefficient 
0-233 in equation (18) is also replaced by 0-204. With all these 
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differences, there is a good deal in hand, and theory predicts that 
oscillations will be maintained under ordinary conditions when all 
the electrons are caught after one excursion. It remains to see how 
subsequent orbits affect the numbers. 


13. Calculation for multiple orbits 

The cosine and sine integrals of §11 are nearly equal when b/a is 
large. If they are put exactly equal, sin(vt,—0@,) alone appears in 
the formulae and much of the work can be done with tables of single 
entry. Orbits which begin with the same value of vt,—6, retain that 
value, and the values of a/a and f, throughout. Only a sketch of the 
calculations is given here. 

We take as before b/a = 500, Y, = 225, and A = 5. Then remem- 
bering that equations (32) to (34) give only amounts acquired in one 
orbit, we have approximately 

cla = 3-148,  1000B/V%, = 22-78, 10008 = 5-498, (40) 
where s is the sum of the sines of the angles vt, —6, for the successive 
orbits up to the beginning of a specified orbit, and c/a, B/W, B are 
the values at the end of that orbit. If the values (40) are placed in 
(24), it becomes a quadratic equation for s in terms of a/a. The two 
values of s are calculated for a/a = 1 to 1-4 at intervals of 0-02, 1-4 
to 2 at intervals of 0-05, and 2 to 12 at intervals of 0-1. Equation 
(25) then gives k* in terms of a/a. In this way a table of the positive 
and negative values of s is constructed for kt = 0 to 1-6 at intervals 
of 0-01. It appears that a/a <1 when 0 < s < 0-00230, so that 
electrons are caught on the filament in that range, the corresponding 
limits of kt being 0 and 0-0180. A table of the filament deflexion 
F = +(180°—2M) is constructed for the same values of k, and 
hence a table of F against s from s = —2-241 through -+-0 to 1-996, 
F going from 29-78° through +180° to —26-32°. This table is not 
at equal intervals of s, but such a table would be inconvenient since 
F varies very quickly when s is small, and more manageable rates 
of change are obtained by taking k* as the variable of uniform 
increase. The formula for deciding whether an electron is caught on 
the anode in any orbit is 


1000(r/b—1) = 5-49s—22-22P cos vt;—22-22Qsinvt,, (41) 


where s stops short at the previous orbit and is zero for the first. 
Tables of 22-22P and 22-22Q for 0, = 0° to 180° at 1° intervals are 
3695.9 Q 
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useful. Taking electrons for which 6, = 0°, 10°,... in the first orbit, 
tables of 6, and r/b-—1 for the subsequent orbits are constructed at 
intervals of 10° in vt,. Capture by the anode is slow, 5 per cent. of 
the electrons describing seven orbits, at which point the calculations 
are stopped. The tables show rapid variation. Few electrons are 
caught on the filament, but many come near it sooner or later, and 
every hear approach means a large deflexion which may defer capture 
by several orbits. No reasonable amount of labour will give the 
current accurately, but a rough estimate is possible. 

Instead of the table in § 12, in which 1000(f,—/,) is tabulated over 
finite ranges, we construct one in which the range is 10° in vt, with 
an assigned middle point; for example, the entry under vt—vt, = 60° 
for any assigned 0, is the difference of the values of f for vt—vt, = | 
and 65°. Consider the electrons emitted within 5° on either side of 


6, = 0 at times within 5° on either side of vt, = 270°. The electron 


6, = 0 has 6, = 55-01° in the second orbit. Its charge is regarded, 
in our rough estimate, as spread over the appropriate range on either 
side of 55-01°, and then the table of 1000(f,—/,) gives the current 
due to the group of electrons at times vt = 0°, 10°,.... The mean 
electron is caught on the anode in the fifth orbit, and its fate is taken 
to decide that of the whole group; that is, the group is considered 
as describing the same number of half orbits as the mean electron. 
The process of calculating the current is straightforward, but pre- 
carious in view of the rapid variations referred to above. From the 
numbers obtained I extract those given on p. 227, in which the first 
line is copied from the table in § 12 (last line), and the others corre- 
spond to orbits after the first. The fundamental term in the con- 
vection current is the real part of 
(0-436— 1-9287) 7 exp wt, 

and since the first term is positive the magnetron in these circum- 
stances will not oscillate. The figures are very different from those 
for a single orbit. Whether the magnetron will maintain a range 
of potentials not including 5 volts can only be decided by detailed 
calculation. 


14. Completion of the theory of an oscillating magnetron 
Let the fundamental term in the convection current be the real 

part of (a,—ib,)/ exp wt, where a, is negative. From § 6 the ampli- 

tude of potential maintained is A = —a,/L/R(K-+Clcosec*vl/c). 
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A value of A is assumed in calculating the current (cf. §13), and 
the last equation gives J, the emission required to maintain that 
potential. A magnetron in which multiple orbits contribute appre- 
ciably to the current is different in principle from one in which they 
do not; for whereas a magnetron with one orbit yields to a com- 
paratively simple linear theory, it is only the effect of the oscillatory 
field on a given orbit which is proportional to A when A/V, is small. 
When the electron approaches the filament it is deflected by an 
amount which varies rapidly with A, and the current in the second 
and subsequent orbits is not proportional to A, except possibly over 
a much smaller range. 

As in § 4 of Note VIII,* the fundamental current need not have 
any particular phase relation with the fundamental voltage. All that 
is allowed for by the slight difference between the frequency of 
oscillation and the natural frequency. FR being small, the harmonics 
are of smaller order. 


15. Strength of oscillations of different frequency in a single- 
anode magnetron 
We return in this section to a single-anode magnetron as simplified 
in Note VIII, taking the undisturbed electrons to describe circular 
orbits and the filament to capture them all on their return. Let the 
radius of the undisturbed cloud be exactly b, the radius of the anode. 
Then if V = Asinvt, equation (3) of Note VIII gives 


Q = Asin(vi—v7/w). 


A gap in the return stream of electrons begins at the anode at time 
a/w, reaches the filament at time 27/w, and lasts for a time z/v. In 


calculating current it is simpler to consider gaps rather than elec- 
trons, so that from equation (5) of Note VIII 
dV log(r,/r.)+log(r3/7,)+... 


‘om rs log(b/a) 


> 


where the gaps extend at time ¢ from r, to 7, 7; to 74, etc. The con- 
vection current is nearly equal to J when r, = a, that is when there 
is a gap at the filament, and nearly zero otherwise. Thus in each 
complete period, say 0 < vt < 27, there is a convection current I 
beginning when vt = 2v7/w or an amount differing from it by a 


* Loc. cit. 213. 
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multiple of 27, and lasting for a half period. If this convection 
current is written 


I s (C,, cos nvt+-D, sin nvt), 
n=0 


70, = —2sin2v7/w and 7D, = 2cos2vz/w. We find, as in §4 of 


Note VIII, that 21 LC cos 2vn|w 
aRK(K+C) ’ 
giving the amplitude of potential maintained in a condenser circuit. 
The condition of maintenance is that cos 2v7/w < 0, so that v/w 
ranges from } to ?, {to j, etc. Taking the factor —cos 2v7/w as the 
measure of the strength of an oscillation, the strength is greatest 
when v/w = }, 3, ete. 


A=— (42) 


16. Some concluding remarks 

We have seen that there is reason for supposing that v = 4w is 
favourable for maintenance in an idealized single anode magnetron. 
This holds for any degree of grazing. It means that the time of the 
slowest oscillation maintained is about four times the time of transit 
of an electron from filament to anode, or twice the time of a complete 
excursion. But the resemblance to mechanical or electrical resonance 
is not very close, and there is no sharp maximum. 

The application of the principle of energy to a problem that has 
already been solved is unnecessary, and I have not considered 
whether the principle of negative work can be improved by allowing 
for the attraction between the electrons and the charge which they 
induce in the plates. But it is instructive to consider the difference 
between low-frequency and high-frequency oscillations. The fre- 
quency may be much less than that at which the wave-length is 
comparable with the distance between the electrodes, and yet it may 
be quite wrong to reckon current by the number of electrons which 
fall on the electrodes. If so, it is idle to inquire on which electrode 
an electron falls, and thus some of the conclusions that are valid for 
low-frequency oscillations cease to hold when the frequency is high. 











A TRANSFORMATION OF EULERIAN 
HYPERGEOMETRIC INTEGRALS 
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1. It was shown by Euler that the hypergeometric series can be 


represented by the definite integrals 


aie Pc) b-1/] __ #)c—b- ree —a 
= pores | (1—t)°-(1—at)-¢ dt, (1) 


0 


F(a, b; c; 2) 


D(C) ‘a— ___ f#)ce—a— __ wt)\— ‘ 
ratece | 1(1—t)°-¢-1(1—2t) dt. — (2) 


0 


Fila, b;¢; 2) = 


These integrals are convergent if, for example, we have |x| < 1 
and c >a,b > 0. If the parameters a, b,c are any real or complex 
numbers such that the integrals are divergent, they can be replaced 
by contour integrals with the same integrands. The paths of inte- 
gration, in the least favourable case, may be double circuits inter- 
lacing the singularities t = 0 and ¢ = 1; and these can sometimes 
be replaced by simple circuits enclosing one or both of these points. 

It was shown by Riemann that a hypergeometric function can 
be represented by forty-eight different Eulerian integrals, which can 
be reduced by an elementary change of the variable of integration 
to one or other of the forms (1) or (2). These cannot in general be 
thus transformed into one another. But Riemann conjectured, and 
Wirtinger subsequently proved, that both forms can be obtained by 
evaluating a certain multiple integral in two different ways. I can 
only mention in passing the striking integral representations of 
ordinary and generalized hypergeometric series recently obtained by 
Erdélyi. Wirtinger’s and Erdélyi’s methods both lie outside the 
scope of this note. 

A few years ago I worked out the analogue of the Riemann- 
Wirtinger theory for the confluent hypergeometric function ,F, (a; c; x). 
This can also be represented by two different classes of Laplace 
integrals, which are derivable from the same double integral. But 
I found incidentally that the circumstances are exceptional when 
Kummer’s confluent hypergeometric equation has a logarithmic 
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singularity; the two types of integrals can then be transformed 
directly into one another by integration by parts. I shall now prove 
that the Eulerian integrals of types (1) and (2), taken along appro- 
priate contours, are also directly transformable into one another, 
when the ordinary hypergeometric equation has a logarithmic singu- 
larity. 

2. If the Riemann equation corresponding to the scheme 


0 oo 1 
Pian B y @ (3) 
a py 
has a logarithmic singularity, we may assume without loss of 
generality that it lies at the point x = 0; for we may permute the 
three singularities by means of the linear transformations 


1 x x—1 
1—z’ x—l’ x 


(4) 


> 


The exponents «a, «’ must differ by an integer; and, since they occur 
symmetrically in the equation, we assume a—a’ = n > 0. We now 
put y = x%(1—a)’y’ to reduce (3) to the standard hypergeometric 


scheme 0 pm 1 


PO a 0 x}, (5) 
l—c 5b c—a—)d | 


where 1—c = a’—a = —n is a negative integer or zero. Thus c is a 
positive integer, while a, b are any real or complex numbers, subject 
only to the condition 
a,b ~ 1,2, 3,...,e—1; (6) 
it is well known that, for these forbidden values, x = 0 is not a 
logarithmic singularity of the hypergeometric equation. The integral 
representation of type (2) is given by a simple contour integral 
(1+ 04 ) 
I= | (1—2t)-%-(1—t)-*-1 dt 
A 

I'(a)T(c—a) 

I'(c) 
whose starting-point A lies in the real interval 0 < ¢ < 1, and where 
we have initially 


am(t) = 0, am(1—t) = 0, jam(1—zt)| < 2. 


= (1—e-?ina) F(a,b;¢;2), (7) 
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If a is an integer not belonging to the forbidden set, the only 
singularity of the integrand is a pole at one or other of the points 
t = 0,1. The path can then be contracted to a point and the integral 
evaluated as a Cauchy residue. The appropriate modifications of 
what follows for this case are left to the reader. 

3. If b is not a forbidden integer (6), we transform (7) by (ec—1) 
successive integrations by parts; the integrated terms cancel around 
the contour, leaving 


r(1 b) (1+,0+) wet 

= — yl-c cs c—b-1 a {2-1 7 c—a-1] J 

I'(c—b) . | (1—2#) qa ( t) | dt. (8) 
A 


But (if a is not a forbidden integer) we have 


oo [e194] = Se Fae $1, 15 15 2) 
die-1  dte-2 a 
I'(a) 
- t@-¢F(a, 1; 1; 
T'(a—c+1) sla 
[(a) 
7 t2-¢( 1 —#)-¢. 
T'(a—c+1) sii . 
Hence vsie uh 


r(i—b)T(a) ae 


— = c—b—lya—c(] __ #)—a 
T'(c—b)l(a—c+1) (1—at)e->-1¢2-€(1 —#) -* dt. 


A 
We now deform the contour into a circle C, |t] = r, where 
l<r<|li/zi, 
and find 
T(1—b)T(a) 


‘= 1 ite—oal 


e-imay1—c | (1—at)e->-l¢e —(t— 1)-¢ dt, 

7 (10) 
where am(t) = 0, am(¢—1) = 0, |am(1—zt)| < 7 at the starting-point 
on the positive real axis. We now write 

e = 1/24, (11) 
which transforms ¢ = 0, 1, 1/x, 00 respectively into s = 0, 1/z, 1, 0. 
The positive circuit C about t= 0 and ¢ = 1 becomes a positive 
circuit about s = co and s = 1/zx, which is topologically the same 
thing as a negative circuit about s = 0 and s = 1, say the circle C in 
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the s-plane described clockwise. Changing the sign and using the 
same circuit as before, we have accordingly 


qi-e [ (1—at)°—>-1¢2-€(¢— 1) -4 dt 


-=f(-J says 


4 { (s—1)¢-13>-1(1 —28)-* ds (12) 


Cc 


(1+,0+) 
— gin(b—c+1) | (1—as)-“s°-1(1—s)¢>— ds, (13) 
A 
where in (12) we take am(s) = 0, am(s—1) = 0, |am(1—zs)| < a at 
the starting-point on the positive real axis, and then deform the 
contour back again into the circuit used in (7). Since we have 
T(1—b) T'(a)e—t74+imb—c+1) 2s P'(a)T(c—a)sin a(c—a) eim(—a4b—c+1) 
T'(c—b)T(a—ce+1) I'(6) (c—b)sin 7b 
_ T@(c—a)(1—e*) 
~~ T(b)T(e—b)(1—e-2**)’ 
because c is an integer, we now get from (7) and (13) the relation 


(1+,0+) 


(1 —xs)—%>-1( 1 —s)eo-1 ds 


(14) 


A 
(1—e-*'7>) T'(b) P(c—b) <i 
a Te) F(a,b;ce;2). (15) 


We have thus derived the representation of type (1) by a direct 
transformation of the type (2). 
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LINEAR PARTIAL DIFFERENTIAL EQUATIONS (I) 
By T. W. CHAUNDY (Ozford) 
[Received 1 April 1938] 

1. An equation of Darboux’s 

Under the term ‘harmonic equations’ Darboux* has discussed a 
class of equations which includes the particular equation 

eV - f{m,(m,+1) ma(m +1) m(ms+1) _m,(m,+1)| V 

oxdy | («+y)? (x—y)? (l—ay)? — (1+ay)? J 


For this equation he gave the solution (symbolic unless m,, m,, ms, 


(1) 


m, are positive integers) 


@\™/ a \m/ a \m/ a \m(t,+t,)2" “i, +4) y(4—ts)) 

ae a ae voll won \ me 7 eat 
(2) 
where t, = 2+y, t, = x—y, t, = 1—xy, t, = 1+2y and ¢, ¢ are 
arbitrary functions. In a formal sense at least this could claim to 
represent the ‘general solution’. My object here is to give the 
fundamental solutiont in terms of which both the general solution and 
the solution of Cauchy’s problem are at once obtainable by well- 

known methods, since the equation (1) is of Riemann’s type 
stap+fpq+yz = 0. 

Such a fundamental solution U(X, Y ; x, y) must have the properties 
that (i) it is a solution of the equation (1) (which is a self-adjoint 
equation), (ii) it reduces to unity along each of the two character- 
istics « = X and y= Y. On this function of the four arguments 














v, y, X, Y impose the transformation 


x = —_&—X)y—Y) an als = ie el 
‘ (@+y(X+Y)’ * (e—y)(X—Y) (3) 
= __(e—X\y—Y) os a —X)y—Y) 
“=~ (1—ay(1— XY)’ 4 > (1-fay)(1+-XY) 
and write 6, = er (r = 1, 2,3, 4). 
ox, 


* Darboux (2), livre 4, chap. 9; in particular 226 (§ 415) equation (47). I 
have modified the notation slightly for my own purposes. 

+ Otherwise the ‘Green’s function’ or the ‘Riemann’s function’. 

t See, for instance, Darboux (2), livre 4, chap. 4 (in particular §§ 358-9) ; 
Goursat (3), chap. 26 (ii), § 498. 
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We then have the equivalence of operators 


@ 1 (X+y 1—Xy Tet 
4(/° 


» 
b+ 9 52+ 84+ 


ax a—X\|a+y 1l—azy 1+ay 
2. ae (r+Y . —#Y 5 1+2¥ . 
ey y—Y\aty ute yet] ny? Thay oT 


In forming the product of ies as we observe that the 
coefficient of 87 is 


(X+y)(z+-Y) . (1), and so on; 
1 


(c—X)(y—Y)(x+y)? (x+y)? 
and that the coefficient of 5, 5, is 
(X+y)(a—Y)+(X—y)(x+Y) 
(x—X)(y—Y)(z?—y?) 
(w—y)(X+Y)+(2-+y)(X—Y) 
(x—X)(y—Y)(x*—y?) 


— — . : , and so on. 


2 


a(x+y)? — x(%—y)? 
There is also a first-order term, namely 
SR, — ie Me 
(ety)? * (e—y)? *" (1—ay)? *  (1+ay)? 
On collecting these results, substituting in (1), and rearranging we 
obtain in virtue of (i), the equation for U 


( + y)? E itd, -38 (5; +8,+83+58,)— i 
v 
~ a pl 8H Be -=8, 2(5,+8,+53+8,)— ma(my+1) [0+ 
+ i = farsa, z 8 s(,+8,+85+5,)—mg(mg+1)|0— 
1 - 
~hcgfetets—2 etn sds+89-mdnct]0 


=Q@, (4) 

We can check the symmetry of this result by noticing that, if we 

change the signs of y, Y in (1), (3), then the suffixes 1, 2, 3, 4 become 

2, 1, 4, 3; if we interchange y, Y with 1/y, 1/Y respectively, the 
suffixes 1, 2, 3, 4 become 4, 3, 2, 1 
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Evidently U will satisfy the equation (4) if it simultaneously 
satisfies the four equations 

§,.(5,+8,+6,+5,)U0 = x,(8,—m,)(6,+m,+1)U (r = 1,2,3,4). (5) 


In virtue of (ii) we look for a common solution of this system which 
reduces to unity when 2,, Xj, %3, x, are all zero. Such a solution 
is given by the multiple hypergeometric series* 


Il (- m,),. (m, +1),,- 


- (6) 
(r3+72+173+74)! 





where, as usual, (¢), denotes ¢(¢+-1)...(¢+-r—1). This series converges 
for sufficiently small x,, i.e. in some neighbourhood of x = X,y = Y. 


2. Some expressions for the function /\ 

This function (6) has some associations with Legendre’s function 
P,. To exhibit its four arguments x, with their associated parameters 
m, I introduce the notationT 


m Mm mM, Mm 
U= A, 1 "a 3 ‘|. 


% Le Bye % 


If an x, say x, (or its associated m) vanishes, the series reduces to the 
corresponding triple hypergeometric series, which I shall denote by 


Mm, Ms, Ms 
Al's | 


a 
similarly defining /\,, and /\,. I shall for brevity write 
[m], = (—m),(m+1),. 


Ale]- SF 


r=0 


Then 


= ,F(—m,m+1;1;2] 


a F,,(1—2z), (7) 
by Murphy’s formula.t 


* Multiple hypergeometric series in two variables (Appell’s hypergeometric 
functions) and their associated partial differential equations are discussed by 
W.N. Bailey (1), chap. 9. 

+ This Chinese character IN may be pronounced (and written) pa. 

t E. T. Whittaker and G. N. Watson (4), § 15.22. 
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a" 
Again A,” “|= p>) 5 nbn b ae" ante 


Now > > = > [---le=o+ > >. 
r=0 8s=9 r=0 r=0 s=1 


The single sum gives P,,(1—2z) as above. The double sum can be 
written 


— [m], [n],a7y® r!(s—1)! 
~ (r!)?s! (s—1)! (r+s)! 


ae > > een [m], [nho'y’ (1 —itye— dt 


(r!)?s!(s—1)! (s—1)! 





0 


- | (Sees ar 


within the region of uniform convergence. Thus 


1 

Al ") = Ra(1—20)—2y { Pg(1—2e +21) Py(1—2yl) dt. (8) 
0 

To put this into a form exhibiting the symmetry we remark that 


P,(1—2a) = [(t—1)P,(1—2a+ 2at)P, (1—2yt) | 


~/@ m | oe 2(1—t)P,, P+ 2y(1—t)P,, P,} dt, 
0 


omitting the arguments of P,,, P, for brevity. Hence 


m? 


1 
Al) = | Pa Py—2e(1—)Pg P,—2ytP, Pd 
0 


i.e. 


1 
Ne 2 ] _ (eZ +u5 +1) | Pa(t— 20+ 20t)P,(1—2yt dt. (9) 
0 


For /\, I use the Dirichlet integral* 


r,!ro! rg! rq! fs 
(Ty +?2+13+14+3)! =|- [LI [ea “ 


* See, for instance, E. T. Whittaker and G. N. Watson (4), 258. 
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where the integration is taken over all non-negative values of the 
t, such that ¢,+¢,+t,+¢, <1. Writing further 


A= 4% +%25— +%35— +27 (11) 
Ox, C9 
we have 
a m, M, Ms, M, 
% fe ty % 
1,!ro!r5!r,! 
r,+r.-t _* 
Ty+TetTst+r,t+3)! 
Substituting from (10) and summing across the integral signs we get 
at length 
Mm, M, Mz, ™ 
1 My My Ms 
Na ae 
Ze ty & 


— (A+1)(A+2)(A+3) ff | | P,,, (1— 2a, ty) P,,(1—2ary ty) X 


X Py (1— 2% ts)P,.(1—2a,t,) dt, dt,dt,dt,, (12) 


ms ms 
where the integral is extended over all non-negative values of ¢, such 
that t+, +t, tt, S a 
As an alternative form to (12) we can write symbolically, 


hao > = TL (1— [ Pi(1—205t,)20,dt,), (13) 


[2% Ly Le XM s=1 


where in the expansion of the product every integral is extended 
over all non-negative values of ¢ such that }¢< 1. This is an 
extension of (8) and is proved similarly. 


3. Some identities 
If in the differential equation (1) we put m, = m, and m, = ms, 
the equation reduces to 


eV (m+ mM: ie 


: 1) 

day Oxy 22— y?)' —x*y?)?’ 

OV  ——- m(m,+1) , m3(ms+1) (14) 
Ca'dy’ = (a’—y’)® * (1—ar'y’)?’ 


9 


‘= y*. This is again of the form (1) with m,, m, zero. 
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Writing similarly X’ = X*, Y’ = Y?, we further find that 
(x*—X*\y*—¥*) _ (@’—X'"\y’—Y') _ ,. 
X,+%,—X,2, = ——,; = SS 7; = I, Bay; 
OES CA ee 
(x*—X*)(y?—Y*) 
~ (l—a®y?)(1— XY?) 
_ @’—Xy'-Y’) 
(l—zy\l—-X) 
Thus the fundamental solution of (14) is 


ae 3h 


% 


%g+%—%s%, = 


la , 
= 7, say. 


Since, as can be proved without difficulty, the fundamental solution 
of an equation of Riemann’s form s+ap-+fq+-yz = 0 is unique, we 
can identify the two fundamental solutions of (14), getting an 
identity of the form* 


ie. aoe el q + 
Ad? Be Xe {|= AS ee Le ice (15) 


Again, if in (14) we put m, = mg, the equation can be written 


2V = m4(my+1)(1—1/2x"*)(1—1/y’?) 


éu'dy’ f(x’ + 1/a')—(y' + 1/y’)P 
eV m,(m,-+ 1) 


(16) 


i.e. = — 
ea" by" (a” 


2? 


—y") 
where x” =2’+1/2’, y” =y'+l1/y’. This again has the form of 
(14) with m, zero. We find that 

(a2” — , X”"\(y”— 7”) 
(x PP sal y")( = x) 


, / , , —s ” 
LetXg—X_X3 = = 2, say. 


Thus further 


m m , m 
° = = Pfs )\i—y)—1}, 17 
by (7). If we combine this with (15) we get for the /\, with four 
equal parameters 


Ale pies; | = Pa(2(1—a)(1—2,)(1—a5)(1—2)—D. (18) 


Z, Be By XM 


* The function IN« is, of course, symmetrical in its four parameters, and 
so (15) applies mutatis mutandis whenever these are equal in pairs. 
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If in these identities we substitute for /\,4, etc., the various 
formulae of §2, we get identities involving integrals of P,,: these I 
have not attempted to verify independently. The simplest of them 
comes by combining (8), (17), whence, writing 1—u, 1—v, $(1+-7) for 
x, y, t, we have 


P,,(2uv—1) 
== P 


m 


1 
(2u—1)—(1—v) | B,fu+(1—u)r}P},fvo—(1—v)r} dr. (19) 
m\ t 


-—% 
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